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We examine the motion in Schwarzschild spacetime of a point particle endowed with a scalar 
charge. The particle produces a retarded scalar field which interacts with the particle and influences 
its motion via the action of a self-force. We assume that the magnitude of the scalar charge is 
small, and that the deviations from geodesic motion produced by the self-force are small. This 
problem is analogous to that of an electric charge moving under the action of its electromagnetic 
self-force, and to that of a small mass moving under the action of its gravitational self-force. We 
exploit the spherical symmetry of the Schwarzschild spacetime and decompose the scalar field in 
spherical-harmonic modes. Although each mode is bounded at the position of the particle, a mode- 
sum evaluation of the self-force requires regularization because the sum does not converge: the 
retarded field is infinite at the position of the particle. The regularization procedure involves the 
computation of regularization parameters, which are obtained from a mode decomposition of the 
Detweiler- Whiting singular field; these are subtracted from the modes of the retarded field, and 
the result is a mode-sum that converges to the actual self-force. We present such a computation 
in this paper. While regularization parameters have been presented before in the literature, there 
are two main aspects of our work that are new. First, we define the regularization parameters as 
scalar quantities by referring them to a tetrad decomposition of the singular field. This is different 
from standard practice, which is to define regularization parameters as vectorial quantities. The 
advantage of dealing with tetrad components is that these, unlike vector components, are naturally 
decomposed in scalar spherical harmonics. Second, we calculate, for any bound orbit around a 
Schwarzschild black hole, four sets of regularization parameters (denoted schematically by A, B, 
C, and D) instead of the usual three (A, B, and C). While only the first three regularization 
parameters are needed to produce a convergent mode-sum, the inclusion of a fourth parameter has 
the practically important consequence of accelerating the convergence. The focus of this paper is 
entirely on the computation of regularization parameters for the scalar self-force. The techniques 
that we introduce in this work are not, however, restricted to this context. They will readily be 
exported to the electromagnetic and gravitational cases, but we leave this generalization for future 
work. As proof of principle that our methods are reliable, we calculate the self-force acting on a 
scalar charge in circular motion around a Schwarzschild black hole, and compare our answers with 
those recorded in the literature. We leave for future work the generalization of this calculation to 
generic orbits. 

I. INTRODUCTION AND SUMMARY 
A. Motivation 

The capture of solar-mass compact objects by massive black holes residing in galactic centers has been identified as 
one of the most promising sources of gravitational waves for the Laser Interferometer Space Antenna [1] . The need 
for accurate templates for signal detection and source identification is currently motivating an intense effort from 
many workers to determine the motion of a relativistic two-body system with a small mass ratio. This is done in a 
treatment that goes beyond the test-mass description in which the small mass moves on a geodesic in the spacetime of 
the large black hole. An additional complication arises from the fact that the treatment cannot rely on a slow-motion 
or weak-field approximation. The accelerated motion of the small mass in the background spacetime of the large 
black hole is governed by the body's gravitational self-force [2, 3], which encodes the influence of the body's own 
gravitational field on its motion. To compute the gravitational self-force acting on a body moving on a bound orbit 
around a Schwarzschild or Kerr black hole is currently the focus of much work; for reviews, see Refs. [4, 5] and a 
special issue of Classical and Quantum Gravity [6] devoted to this topic. 

The complexities associated with the computation and interpretation of the gravitational self-force have motivated 
the formulation of educational toy problems. These have the advantage of being much simpler to deal with, but they 
should nevertheless capture the essential physics of self-forced motions in curved spacetime. One such toy problem, 
which is based on real- world physics that is interesting in its own right, is the motion of an electrically charged particle 
in curved spacetime, subjected to a self-force produced by its own electromagnetic field. The foundations for this 
problem were laid more than 45 years ago by DeWitt and Brehme [7] (their work was later corrected by Hobbs [8]). 

Another such problem is the motion of a particle endowed with a scalar charge; here the foundations were provided 
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more recently by Quinn [9] . In spite of its academic nature (there are no known macroscopic particles endowed with 
a scalar charge), the scalar self- force problem is a useful toy problem because of its relative technical simplicity, and 
because the motion of a particle under the influence of its scalar self- force is expected to resemble closely the motions 
obtained in the other contexts (electromagnetic and gravitational). The self- forced motion of a scalar charge in 
Schwarzschild spacetime is the focus of the work presented in this paper. We exploit this simple situation to introduce 
computational techniques that are required for the concrete evaluation of the self-force. These techniques, however, 
are not limited to the context of the scalar self-force, and they will readily be exported to the electromagnetic and 
gravitational cases. We leave this generalization for future work. 



B. The problem 

Consider a particle of mass m and scalar charge q moving on a world line 7 of the Schwarzschild spacetime. The 
world line is described by the parametric relations z a (r) in which r is proper time. The particle produces a retarded 
scalar potential <I> that satisfies the inhomogeneous wave equation [9] 

n$(x) = -4nq ( 5 4 (x,z)dT, (1.1) 
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in which □ := (/"^V^V^ is the covariant wave operator, and 8i{x, z) is a scalarized Dirac distribution with support 
on the world line. Except for the factor — A-k inserted for convenience, the right-hand side of Eq. (1.1) is the particle's 
scalar charge density. 

The retarded potential <&{x) produces a field ^ a ( x ) '■— V Q $(a;) that acts on the particle and influences its motion. 
As shown by Quinn [9], the particle's acceleration is proportional to the components of $ Q that are orthogonal to the 
particle's velocity vector u a := dz a /dr, 



ma 



q(g^ + u a u^p(z), (1.2) 



where a a := Du a /dr is the particle's acceleration vector, the covariant derivative of the velocity vector along the 
world line. Quinn also showed that the longitudinal component of the retarded field is responsible for a change in the 
particle's inertial mass, 

nun 

?P = -qu a $ a (z). (1.3) 
dr 

This effect was explored in cosmological situations in Refs. [10, 11]. 

Equations (1.2) and (1.3) have only formal validity because the singular field $ Q (x) diverges as x — > z: the field of 
a point charge is necessarily infinite at the position of the particle. Quinn [9] was able to regularize these equations so 
as to produce meaningful equations of motion for the charged particle. In breakthrough work that plays a central role 
in this paper, Dctwcilcr and Whiting [12] showed that Quinn's regularization procedure amounts to a decomposition 
of the retarded potential into uniquely defined singular and regular potentials, 

*(a;) = ® s {x) + <S> R (x), (1.4) 

with $ s denoting the singular potential and <& R the regular potential. As was shown by Dctwcilcr and Whiting, the 
singular potential possesses the following properties: (i) it satisfies the same wave equation as the retarded potential, 
that is, it is a solution to Eq. (1.1); (ii) it displays the same singularity structure as the retarded potential near the 
particle's world line; and (iii) it does not exert a force on the point charge. These properties make the decomposition 
of Eq. (1.4) unique. The regular potential, on the other hand, possesses the following properties: (i) it satisfies a 
homogeneous version of Eq. (1.1), with a zero right-hand side; (ii) it is smooth on and near the particle's world line; and 
(iii) it alone determines the self-force acting on the particle. The conclusion, therefore, is that the actual equations of 
motion for the particle are Eqs. (1.2) and (1.3) with the regular field $ R := V Q $ R substituted in place of the retarded 
field <!> Q ; this conclusion is in full agreement with Quinn's earlier work. The Detweiler- Whiting decomposition also 
plays an essential role in the electromagnetic and gravitational self- force problems. 

In Schwarzschild spacetime it is computationally advantageous to solve Eq. (1.1) after decomposing $(x) in spherical 
harmonics Yi m (6, (f>). Adopting the usual Schwarzschild coordinates [t, r, 8, <f>], we would write 

=2$j m (t,r)y Im (M) (1-5) 

lm 
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and call &i m (t, r) the spherical-harmonic modes of the retarded potential; the right-hand side involves a sum over 
all integers I = 0, 1, 2, ■ ■ ■ , oo and a nested sum over all integers m = —I, — I + 1, • • • , I — 1, 1. After performing the 
decomposition of Eq. (1.5), Eq. (1.1) turns into a two-dimensional wave equation for each mode function $( TO (t,r); 
this equation is displayed in Eq. (1.12), below. 

The spherical-harmonic modes of the retarded potential give rise to quantities & a i(x) := V Q J2 m ^Im^mi such. 
that the retarded field can be expressed as the mode-sum 



* a (x) =]r$ Q ^)- (1.6) 

i 

We call these quantities the multipole coefficients of the retarded field. Each $ Q ;(x) is bounded in the limit x — > z, 
in spite of the fact that the retarded field is infinite on the particle's world line. (The multipole coefficients are 
discontinuous at x = z.) The mode-sum of Eq. (1.6), of course, does not converge when the retarded field is evaluated 
on the world line. 

The failure of Eq. (1.6) to converge when x — z is exactly compensated for by the failure of 

5>a«(s) (1.7) 
I 

to converge, because (as was noted previously) the retarded and singular fields share the same singularity structure 
near the world line. Here, := V Q ^ m ^fmXlm with $f m denoting the spherical-harmonic modes of the singular 
potential. The regular field can thus be expressed as 

**(z)=\hnY^\**lW-*ll(*j\> ( L8 ) 

x—>z ' L 
I 

in terms of a converging mode-sum. The limiting procedure involved in Eq. (1.8) is introduced to handle the (shared) 
discontinuity of the multipole coefficients <J> Q /(x) and ^ t (x) at x = z. After subtraction the multipole coefficients of 
the regular field are smooth at x = z, the mode-sum converges, and the limit can be taken without difficulty. 

The prescription contained in Eq. (1.8) becomes a practical method to evaluate <I>^ (z) — and therefore the self-force 
acting on the scalar charge — when one can actually compute the quantities ; (x) for a field point x close to the 
world line. As we shall review below, this computation is possible because the singular field 3>|(a;) is known in a 
neighborhood of the world line; it can be expressed as a Laurent expansion in powers of the distance to the world 
line. The end result for the multipole coefficients takes the schematic form 



{l + \)A a + B a ' C 



(1.9) 



in which the quantities A a , B a , and C a , known as regularization parameters, are independent of I but depend on 
the state of motion of the particle at z. Notice that the sum over the A a term would diverge quadratically, that 
the sum over the B a term would diverge linearly, and that the sum over the C a term would diverge logarithmically. 
The remaining terms, those designated with (•••), lead to a converging sum that evaluates to zero by virtue of the 
Detweiler- Whiting axiom, according to which the singular field docs not produce a force on the particle. This ensures 
that after removal of the A a , B ai and C a terms, the mode-sum of Eq. (1.8) will converge to the correct value for 

Our central task in this paper is the computation of regularization parameters for a scalar charge moving on a 
bound orbit around a Schwarzschild black hole. We leave for future work the completion of a calculation of the actual 
self-force. This would involve, over and above the work presented here, a numerical determination of the retarded 
field produced by a scalar charge moving in Schwarzschild spacetime. 

We consider the particle's charge q to be small, and we have in mind a perturbative implementation of Eqs. (1.2) and 
(1.3). In a zeroth-order approximation, the particle is taken to have a constant inertial mass and to move on a geodesic 
of the Schwarzschild spacetime. In a first-order approximation, the regular field is computed for this geodesic 
motion and substituted on the right-hand sides of Eqs. (1.2) and (1.3). This iterative process could be continued, 
but we suppose that q is sufficiently small that the process can be stopped after a single iteration. In this regime 
the self-force can be computed while assuming that the motion is geodesic. We adopt this small-charge approximation 
here, noting that it reflects the spirit of the small-mass-ratio approximation in the gravitational problem. The 
assumption leads to much simplification; the computation of regularization parameters for accelerated particles would 
be significantly more involved. 
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C. 



Past work 



We are not the first researchers to define and compute regularization parameters for the mode-sum computation of 
self-forces in curved spacetime. In fact, the literature is vast and the field has already reached a fairly mature state. 
But as we argue below, we believe that this work (and our promise for extensions toward the electromagnetic and 
gravitational problems) is a significant addition. 

The main ideas behind the mode-sum regularization of self-forces were first introduced by Barack and Ori in 
their pioneering work [13-15], which was later perfected and extended by an Israeli- Japanese consortium including 
Barack, Ori, Mino, Nakano, and Sasaki [16-19]. These authors computed regularization parameters for the mode-sum 
evaluation of the (scalar, electromagnetic, and gravitational) self-force acting on a particle moving on any bound 
geodesic of the Schwarzschild spacetime. Barack and Ori [20] were then able to extend these results to the Kerr 
spacetime. This early work on regularization parameters is nicely summarized in Ref. [21]. 

The (so-called early) work reviewed in the preceding paragraph was carried out before the discovery by Dctweiler 
and Whiting of the retarded field's decomposition into uniquely identified singular and regular pieces [12]. It relied 
on an alternative decomposition, in terms of "direct" and "tail" pieces, which did not come with the same degree 
of mathematical elegance. For example, unlike the regular field which is smooth everywhere in a neighborhood of 
the world line, the tail part of the retarded field is not diffcrcntiablc on the world line. The Detweiler- Whiting 
decomposition provides a much sounder foundation for the definition and computation of regularization parameters. 
In follow-up works, Detweiler, Messaritaki, and Whiting [22] carried out such a computation for a scalar charge moving 
on a circular orbit around a Schwarzschild black hole, and Kim [23] extended these results to a generic orbit. 

A number of works present explicit computations of the self-force (for various charges undergoing various motions 
in various spacetimes) by mode-sum techniques. Burko computed the self-force acting on an electric charge in circular 
motion in Minkowski spacetime [24] . He and his coworkers also considered scalar and electric charges kept stationary 
in a Schwarzschild spacetime [25], in a spacetime that contains a spherical matter shell (Burko, Liu, and Soen 
[26]), and in a Kerr spacetime (Burko and Liu [27]). In addition, Burko computed the scalar self- force acting on a 
particle in circular motion around a Schwarzschild black hole [28]. This calculation was since revisited by Detweiler, 
Messaritaki, and Whiting [22], as well as Diaz-Rivera, Messaritaki, and Whiting [29], who also considered the case 
of slightly eccentric motion. Barack and Burko dealt with a particle falling radially into a Schwarzschild black hole, 
and evaluated the scalar self- force acting on this particle [30]; Lousto [31], and Barack and Lousto [32], computed the 
gravitational self-force for radial infall. 



In this paper we present a new calculation of regularization parameters for the mode-sum evaluation of the scalar 
self-force acting on a particle moving on a bound geodesic of the Schwarzschild spacetime. Our calculation is new in 
two main respects. 

First, we define the regularization parameters as scalar quantities by referring them to a tetrad decomposition of the 
singular field. In contrast, the original parameters of Eq. (1.9) are vectorial quantities that refer to <f>^, the vectorial 
components of the singular field. The idea here is to introduce a basis of orthonormal vectors e?^ (x) at every point of 
the Schwarzschild spacetime; the superscript a is the usual vectorial index, and the subscript (fi) = {(0), (1), (2), (3)} is 
a label that designates an individual member of the tetrad. The vectors satisfy <?a/3e^e^ = = diag[— 1, 1, 1, 1]. 

The four quantities ^ := ef^a are the frame components of the singular field, and these are scalar functions of 
the spacetime coordinates. 

The advantage of introducing the tetrad and the associated decomposition of vector fields is that each frame com- 
ponent •. (x) is a scalar function that can naturally be expanded in scalar spherical harmonics. This is quite unlike 

the vector $>^(x), which could be expanded elegantly in vectorial harmonics (a procedure that has not been adopted 
in the self-force literature) or inelegantly in scalar harmonics (as was done in all previous works on regularization 
parameters). By introducing the tetrad we are able, in this work, to provide an elegant definition for the regularization 
parameters. Our specific choice of tetrad will be specified below. 

Second, we calculate an additional set of regularization parameters in order to accelerate the numerical convergence 
of the mode-sum. Together with our implementation of the tetrad decomposition, this amounts to replacing Eq. (1.8) 



D. This work 



by 




(1.10) 
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and Eq. (1.9) by 

r n,..-, n,.. i 

(1.11) 



(/ + 2) A (li) + B M 



(j+i) (Z- 3) 



The regularization parameters A^, B^, and C^) have already appeared (in vectorial form) in Eq. (1.9); the 
regularization parameters are new. We calculate all of these for a scalar charge moving on a bound geodesic of 
the Schwarzschild spacetime. 

We note that the vectorial parameters D a were computed by Detweiler, Messaritaki, and Whiting for the special 
case of circular motion in Schwarzschild spacetime [22], and that they were computed by Kim for generic orbits; Kim's 
results are recorded in his PhD dissertation [33], but they have yet to appear in the peer-reviewed literature. We note 
also that because the vectors e"^ contain an angular dependence, the operations of [multiplication by a tetrad vector] 
and [extraction of multipole coefficients] do not commute: (e"^^ a )l ^ e^)(^a);; as a consequence, our expressions 
for -B( M ), C(fi)-, and D(n) cannot be compared directly to those for A a , B a , C a , and D a that have appeared in 

the literature. 

It is important to mention that the sum over / of the D^) term in Eq. (1.11) is actually zero, because Ym^oK^ ~ 
|)G + = 0- And as we have seen, the same statement applies to the remaining terms in Eq. (1.11), those 

designated by (•••). These terms, therefore, do not contribute to the final value of <&?n(z) when the sum is evaluated 
in full, from / = to I = oo. Nevertheless, the term does play a useful role when the sum is truncated to 

some finite upper bound I = l max : It produces a significant acceleration of the sum's convergence. This property has 
practical importance, because to truncate the sum over I is a computational necessity. 

We compute the regularization parameters A^, B^, C( M ), and by importing many techniques from the 

literature. In fact, an advantage of arriving late into this business is that we can pick and choose, from a variety 
of sources, the computational methods that are the most compelling. Thus, for example, in Sec. II we develop a 
covariant local expansion for 3>„(a;) that is inspired by Mino, Nakano, and Sasaki [19], but which rests on the more 
secure Detweiler- Whiting decomposition of the retarded field into singular and regular pieces [12]. As another example, 
in Sec. V and the Appendix we employ the rotated angular coordinates (a,/3) of Barack and Ori [17], but we do so 
without having to deal with vector components that are singular at the point of evaluation (a = 0). As a third 
example, in Sec. V we routinely make substitutions such asa-> y/2 — 2 cos a + 0(a 3 ) to turn a function of the angles 
that is well-defined only in a neighborhood of a — to another function that is well-defined on the entire sphere. We 
stole this powerful idea from Detweiler, Messaritaki, and Whiting [22], but we implement it without making contact 
with the Thorne-Hartle-Zhang coordinates [34, 35]; these form an important part of their analysis, but they play no 
role here. 



E. Prescription 

We summarize our main results in the form of a detailed prescription for the concrete evaluation of the scalar 
self-force. We recall that the force is acting on a particle moving on a generic orbit around a Schwarzschild black hole, 
and that it is evaluated within the small-charge approximation described near the end of Sec. I B. 

First step: Integrate the wave equation. The first task that must be accomplished is to solve Eq. (1.1) for 
the retarded potential $(x). This is best accomplished by decomposing the potential in spherical harmonics, as in 
Eq. (1.5). Each mode <&; m (t,r) of the retarded potential satisfies the reduced wave equation 



dt 2 drl 



1(1 + 1) 2M 



Y^ lm ) = -^Y lm (%,0)c- im ^ 5(r - r(t)), (1.12) 



in which the right-hand side is the reduction of the scalar charge density — Anq J 8a(x, z) dr. This equation must be 
solved while imposing ingoing-wave boundary conditions at the black-hole event horizon, and outgoing- wave boundary 
conditions at infinity. The modes are complex functions, related to each other by = (— l) m $; m , in which an 

overbar indicates complex conjugation; this condition ensures that the scalar potential of Eq. (1.5) is real. 

In Eq. (1.12), the reduced wave operator is written in terms of the tortoise coordinate r* = J / _1 dr = r + 
2Mln(r/2M — 1), where M is the black-hole mass and 

2M 

/:=!-— (1-13) 



6 



The reduced charge density depends on the functions r(t) and <p(t), which give the coordinate positions of the world 
line. These are obtained by solving the geodesic equations, 



t 



E 



1 - 2M/r' 



1+^7- 



L 

72' 



(1.14) 
(1.15) 
(1.16) 



in which an overdot indicates differentiation with respect to proper time r. The constant E is the particle's conserved 
energy per unit rest-mass, and the constant L is the conserved angular momentum per unit rest-mass. It is assumed 
that the particle is moving with 9 = ^ in the black hole's equatorial plane, and the right-hand side of Eq. (1.12) also 
involves w* — i = E/f, the time component of the particle's velocity vector. 

Equations (1-12) and (1.14)-(1.16) can be integrated with any reliable numerical method, and this procedure returns 
^im{t, r) for selected values of I and m. From this we extract the quantities 



d d 
$lm{t,r + ), —$i m (t,r + ), — $j m (i,r+), 



(1.17) 



which are evaluated at r 



r(t) + A, slightly away from the radial position of the particle at time t. The symbol 



A represents a small radial displacement, which can be either positive or negative. While the first two quantities 
listed above are actually continuous across r = r(t), the third quantity is discontinuous, and by evaluating it at r = r+ 
we make its value unambiguous. 

Second step: Convert the modes. The first step provides us with the spherical-harmonic modes of the retarded 
potential The computation of the self-force, however, as implemented in Eq. (1.10), requires the spherical-harmonic 
modes of := e? )V a <&, the frame components of the retarded field in the selected basis of orthonormal vectors. 
Our choice of tetrad is motivated in Sec. IV; in the usual ordering [t, r, 9, <p] of the Schwarzschild coordinates we have 



e (o) 



6 (2) = 



0,0,0 



sm< 



0, \/J sin 9 cos <f>, - cos 9 cos v , 

r r sin 



0, \ff sin sin <f>, — cos 9 sin <j), ' 



r sin 9 



"(3) 



(1.18) 
(1.19) 
(1.20) 
(1.21) 



0, \ff cos 9, sin 9, 

r 

In practice it is useful to introduce, as substitutes for e™^ and e" 2 ^, the complex combinations e" ± ^ := e"-^ ± ie" 2 ^, or 

(1.22) 



e (±) - 



i— 1 +io ±i< ^ 

0, Jf sin Oe^*, - cos 0e ±ifl \ - 

r r sin 9 



As shown in Sec. IV, the spherical-harmonic modes < I ) ( M )( m (i, r) are given in terms of $; TO (t, r) by 



(0)lr, 



(1.23) 



(+)Zm 



'(Z + m— 1)(Z + m) 
(2Z — 1)(2Z + 1) 



or r 



+ 



'(Z-TO + l)(Z-TO + 2) 

(2/ + 1)(2Z + 3) 



Z + 2 



77# + 

or r 



i+l,m-l) 



(1.24) 



(-)lm — 



! (l-m-l)(l-m) 
(2Z — 1)(2Z + 1) 



Z - 1 



'(Z + m+l)(Z + m + 2) 
(2Z + 1)(2Z + 3) 



V7 



r 

d_ 

dr 



$J-l,m+l 

Z + 2 S 



i+l,m+l) 



(1.25) 
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(3)ir, 



+ 



' (/-m)(Z + ro) / ,-0 
(2/-l)(2i + l) ^ V/ 9r 

' (I - m+ + m + 1) 
(2Z + l)(2Z + 3) 



Z-l 



i-i,' 



/7 9 , * + 2 \ii 



(1.26) 



These can be evaluated at r = r + by substituting the values extracted in Eq. (1.17). 

Third step: Regularize the mode sum. The spherical-harmonic modes $( M )i m (f, r) give rise to the multipole 
coefficients of the retarded field, which are defined by 



*(„)i(t,r,0,0) := ^ $ M ; m (i,r)y; m (6»,^). 



(1.27) 



The frame components of the retarded field are then given by 



(1.28) 



This sum can be evaluated at r = r+ := r(t) + A, 6 = §, and = (£>(£). The sum would diverge in the limit A — > 0, 
but it is to be regularized as in Eq. (1.10), which we rewrite as 



( /+i) (/-!)(/ + §) 



+ 



after also involving Eq. (1.11). 

We compute the regularization parameters in Sec. V of this paper. We find 



(1.29) 



1(0) 



A 



(+) 



Vf (r 2 + L 2 ) 



sign(A), 



Vf (r 2 + i 2 ) 



sign(A), 



A {3 ) = 0, 



(1.30) 

(1.31) 
(1.32) 



where f := 1 — 2M/r and sign(A) is equal to +1 if A > and to —1 if A < 0. We have, in addition, ^4(_) = 
A {1) = Re[A (+) ], and A {2) = lm[A {+) }. 



We also find 



B, 



Err 



(o) 



Err 



Vf (r 2 + L 2 ) 3 / 2 2Vf (r 2 + L 2 ) 3 / 2 ' 
B (+) = e^flfo-iBfo), 



(+) 



Vf(r 2 +L 2 ) 3 / 2 2rVr 2 + L 2 

Vf-i 



+ 



B l+) - - 



2Vf{r 2 + L 2 f/ 2 rVr 2 + L 2 . 
(2-Vf)r . , (2-Vf)r 



+ ■ 



2LVr 2 + i 2 Vf 2iVr 2 + L 2 V~f 



B, 



(3) 



0. 



(1.33) 
(1.34) 

(1.35) 

(1.36) 
(1.37) 



We have, in addition, £>(_) = £(+), £>(i) = Re[B(+)] = costp + sin 93, and B( 2 ) = Im[i?( + )] = B^ simp — 

BJ + s COS if. 

We have introduced the (rescaled) elliptic integrals 



S 



:= - f ' (l-/csin 2 ?A) 1 / 2 rfV = F(-i,i;l;fc) 

7f Jo 



(1.38) 



8 



and 



2 r /2 

X ■=- \ (i-fcsin 2 vO~ 1/2 # = -F(ii;i;£Oi 

n Jo 



(1.39) 



in which k := L 2 /(r 2 + L 2 ). As indicated in Eqs. (1.38) and (1.39), the elliptic integrals can also be expressed in terms 
of hypergeometric functions. 
We also find 



and 



D 



(0) 



D (+) = 



C W = o 



£r 3 (r 2 - L 2 )r 3 E(r 7 + 30Mr 6 - 7L 2 r 5 + 114ML 2 r 4 + 104MLV + 36ML 6 )r 



(1.40) 



2y/f (r 2 + L 2 ) 7 / 2 



m A ^(r 2 + L 2 fl 2 



+ 



£r 3 (5r 2 - 3L 2 )r 3 £(r 5 + 16Mr 4 - 3L 2 r 3 + 42AlX 2 r 2 + 18AlX 4 )r 



16^ (r 2 + L 2 ) 7 / 2 



r 3 (r 2 -L 2 )r 4 



(+) 



16r 2 V?(r 2 + L 2 ) 5 / 2 
(3r 7 + 6Mr 6 - L 2 r 5 + 31ML 2 r 4 + 26A4X 4 r 2 + 9A4X 6 )r 2 



(1.41) 
(1.42) 



+ 



2Vf(r 2 + i 2 ) 7 /2 4(r 2 + i 2 ) 3 / 2 4r 4 Vf (r 2 + L 2 ) 5 / 2 

(3r 7 + 8Mr 6 + L 2 r 5 + 26ML 2 r 4 + 22ML 4 r 2 + 8ML 6 )Vf r 3 + 2Afr 2 + 4ML 2 



16r 6 (r 2 + L 2 ) 3 / 2 



+ 



r 3 (5r 2 -3L 2 )r 4 



+ ■ 



8r 4 Vr 2 + L 2 
(7r 5 + 12Mr 4 - L 2 r 3 + 46ML 2 r 2 + 18AfL 4 )r 2 



+ 



16r 2 V?(r 2 +L 2 )5/2 



D 



(+) 



(3) 



+ 



+ 



0. 



Wf (r 2 + L 2 ) 7 / 2 8(r 2 +L 2 ) 3 / 2 
(7r 5 + 6Mr 4 + 6L 2 r 3 + 12A4X 2 r 2 + 4ML 4 )v / f 

16r 4 (r 2 +L 2 ) 3 / 2 ' SrVr^+T 2 ". 

r 2 (r 2 - 7L 2 )(v / f- 2)r 3 (2r 7 + Afr 6 + 5L 2 r 5 + 10ML 2 r 4 + 29ML 4 r 2 + UML 6 )r 
16LVf(r 2 + L 2 fl 2 8r 5 L(r 2 +L 2 ) 3 / 2 
(r 5 - Mr 4 + 4L 2 r 3 - 5ML 2 r 2 + 2ML 4 )r] 



(1.43) 



4r 3 iV? (r 2 + L 2 ) 3 / 2 
r 2 (r 2 - 3L 2 )(V? - 2)r 3 (4r 5 + 2A/r 4 + 7L 2 r 3 + 10A1XV + 14ML 4 )r 



+ 



16LVF (r 2 + L 2 ) 5 / 2 
(2r 3 - 2Mr 2 + 5L 2 r - 8ML 2 )r 



16r 3 L(r 2 + L 2 ) 3 / 2 



SriVf (r 2 + L 2 ) 3 / 2 



(1.44) 
(1.45) 



We have, in addition, £>(_) = -D(i) = R c [-D(+)] = cos< P + sin<p, and D( 2 ) = Im[Z?( + )] = £>( + ) sin 95 — 

cos <p. 

Fourth step: Construct the self-force. The mode-sum of Eq. (1.29) can now be evaluated. Thanks to 
the presence of the regularization parameters D^, the sum converges quickly to a precise estimate for the frame 
components $j^(z). The vector field := V Q $ R is related to these by 



-*(0)e(0)a + 2 $ f+) e (-)" + 2*f-) e (+)« + ^) e (3W 



After involving Eqs. (1.18)-(1.22) and evaluating this at r = r(i), 6 = ^ , and <p — <p{t), we obtain 

<3> R = i f <f> R ,a — ' v -I- .p' 1 ''! 

o,/fV (+) (-) /' 



2Vf 

-r*fs). 

-^$f +) e^-$^ 



(1.46) 

(1.47) 
(1.48) 
(1.49) 
(1.50) 



These, finally, can be substituted into Eqs. (1.2) and (1.3) for a concrete evaluation of the scalar self-force. 
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A 



•••• 



•••••• 



FIG. 1: Plot of several regularized versions of |Re$( + );| as a function of multipole order I. The different curves are described 
in the text. 



F. Case study: Particle on a circular orbit. 

The prescription detailed in the preceding subsection will be fully implemented in a future publication. To convince 
ourselves that it actually works, we carried out the computations for the special case of a scalar charge moving 
on a circular orbit around a Schwarzschild black hole. Our results are not new: They reproduce some already 
obtained by Burko [28], Detweiler, Messaritaki, and Whiting [22], as well as Diaz-Rivera, Messaritaki, and Whiting 
[29]. Nevertheless, we present them here (without derivations) because they constitute a proof of principle that the 
prescription is valid. 

We place a scalar charge on a circular, geodesic orbit at a radius ro = 6M (this is the innermost stable circular 
orbit). We go through all the steps listed in Sec. I E and compute ^Q^y the frame components of the regular field 

evaluated at r = r(t) := r , 9 = ^, and 4> = : = Qt, where = M/r\ is the particle's angular velocity. Without 
loss of generality we evaluate the scalar field $^ at t = 0, so that ip = 0. This is related to the frame components by 

* t R = V7b*fo). $? = -^Re$f +) , $* = 0, $«=r Im$f +) , (1.51) 

where / = 1 — 2M/r . The regular field satisfies the helical condition 

$R + Q$R = 0, (1.52) 

and we find that the real and imaginary parts of fully determine the scalar self-force. 

In Fig. 1 we plot several regularized versions of |Re$( + ) J |, the absolute value of the real part of <&(+);, the multipole 
coefficients of the frame component $( + ) of the retarded field; we plot this on a logarithmic scale as a function of I, in 
the interval < I < 40. The upper curve (in open triangles) represents the unrcgularizcd multipole coefficients of the 
retarded field; we sec that this function grows linearly with I, so that the sum of its terms diverges. The second curve 
(in open squares) is what is obtained after subtracting (I + ^)Ai + \ from the first curve; we see here that the curve is 
approximately constant, so that its sum also diverges. The third curve (in open diamonds) is what is obtained after 
subtracting Bi + \ and Cr + \/(l + \) = from the second curve; this produces a function that decays as 1/Z 2 , and this 
leads to a converging sum. The convergence is accelerated, however, with the fourth curve (in solid circles), which is 
obtained after subtracting £>(+)/[(/ — + |)] from the third curve; this produces a function that decays as l/l 4 . 

Figure 1 constitutes a very robust test of our numerical and analytical computations: Any error would give rise to 
a gross violation of the properties listed above. For example, an error in the analytical form of would produce a 
curve in full circles that would still decay as l/l 2 instead of the observed l/l 4 . Similarly, a coding error would return 
a retarded field whose singularity structure would not be compatible with the rcgularization parameters, and this 
would again produce very visible effects in Fig. 1. We have tested this observation by deliberately inserting errors at 
various places in our numerical code. The results follow expectations and convince us that the prescription is valid. 
(We have not yet tested the r ^ sector of the prescription.) 

In Fig. 2 we plot |Im$( + ) ; |, the absolute value of the imaginary part of *&(+)/; again we plot this on a logarithmic 
scale as a function of I, in the interval < / < 40. This curve requires no regularization: the imaginary parts of ^4(+), 
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FIG. 2: Plot of |Im$( + );| as a function of multipole order I. 



and £>(+) all vanish when r = 0. (Recall that we evaluate the field at t = and tp := fit = 0.) Here we see 
that |Im$( + );| decays exponentially as a function of I (approximately as e~ 1 / 2 ) until round-off errors start to dominate 
when I is approximately equal to 20. This curve, of course, produces a rapidly converging sum. 
Our final numerical results are 

M 2 

— ®? ^ 3.60907254 x 10" 4 , (1.53) 
M 2 

— ~ 1.67730 x 10~ 4 , (1.54) 

~ -5.30423170 x 10" 3 , (1.55) 

with the number of significant digits reflecting our best estimation of the code's numerical accuracy (the last digit 
is uncertain). The least accurate number is for which is obtained after several rounds of regularization. Our 
numbers are consistent with results obtained by Diaz- Rivera, Messaritaki, and Whiting [29]: In their Table I they list 
(M 2 /q)<Pf = 1.6772834 x 10~ 4 for r = 6M. 



G. Organization of this paper 



The chain of calculations that lead to the prescription detailed in Sec. I E is a long one, and it occupies the remaining 
sections of the paper. Here is how the rest of the paper is organized. 

We begin in Sec. II with the development of a covariant local expansion of the singular field <I>^ (x) in the vicinity 
of the particle's world line. The expansion is based on the assumption that the scalar charge follows a geodesic 
of a vacuum spacetime, but it is otherwise general; we do not yet, at this stage, assume that the metric is given 
by the Schwarzschild solution. The expansion must be carried out to a sufficient degree of accuracy to permit the 
determination of all four regularization parameters. Such an accurate expansion has never appeared in the literature, 
and we present it here for the first time. 

In Sec. Ill we convert the covariant expansion into an explicit coordinate expansion that can be evaluated for 
any spacetime whose metric is expressed in any coordinate system. The methods by which we obtain the covariant 
and coordinate expansions rely heavily on the general theory of bitensors. These were introduced by Synge [36] and 
DeWitt and Brehme [7], and the theory is conveniently summarized in Poisson's contribution to Living Reviews in 
Relativity [4]. This last paper is an essential resource for the calculations presented in Sees. II and III, and we will 
repeatedly refer to it as LRR. 

In Sec. IV we motivate our choice of tetrad e"^. We also work out the relationships listed in Eqs. (1.23)-(1.26) 
between <&(^); TO , the spherical-harmonic modes of ( E>( M ) := e? >Va$, and $; m , the modes of the scalar potential. As 
we explain in this section, the tetrad is selected so as to produce a simple relationship in which $( /i )/ m is linked to 
the neighboring modes $z±i iT „ and $i±i im ±i only; other tetrads would lead to more complicated couplings and are 
best avoided. 
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In Sec. V we compute the regularization parameters Ar^\, Br^, CV^) = 0, and D r^y We do this by expanding 
our local coordinate expansion for <f>^ (a;) in Legendre polynomials and showing that the result takes the form of 

Eq. (1.11). To perform the Legendre decompositions we rely on techniques imported from Dctwciler, Messaritaki, and 
Whiting [22] ; these are summarized in the Appendix. 

Many calculations presented in this paper are extremely tedious and could not have been carried out with pen and 
paper. We relied heavily on the symbolic manipulator GRTensorII [37] working under Maple. Throughout the 
paper we use geometrized units in which G = c = 1, and we adhere to the sign conventions of Misncr, Thorne, and 
Wheeler [38]. 



II. COVARIANT LOCAL EXPANSION OF THE SINGULAR FIELD 



A. Singular field 

The singular part <3? s of the retarded potential $ produced by a point scalar charge q moving on an arbitrary 
world line of an arbitrary curved spacetime was first correctly identified by Detweiler and Whiting [12]. As they have 
shown, the singular potential possesses the following properties: (i) it satisfies the same wave equation as the retarded 
potential, Eq. (1.1); (ii) it displays the same singularity structure as the retarded potential near the particle's world 
line; and (iii) it does not exert a force on the point charge. The scalar self-force acting on the charge therefore results 
from the sole action of the regular potential $ R = $ — $ s , which is smooth on the world line. The self- force is 
proportional to $ R := V a $ 11 , and it can be calculated by first computing $ a := V a $, then removing from this the 
singular part := V Q $ S , and finally evaluating the result at the position of the particle. 

Our task in this section is to develop a covariant expansion of the singular field <j>^ in powers of e, a book-keeping 
quantity that loosely represents the distance between the field point and the world line. (The distance to the world 
line will be defined precisely below.) As was justified near the end of Sec. I B, we shall restrict our attention to the 
case of a scalar charge q that moves on a geodesic of a vacuum spacetime; the charge's acceleration vector and the 
spacetime's Ricci tensor will therefore be set equal to zero. The expansion begins with a term of order e~ 2 and we 
shall keep it accurate through order e, neglecting terms of order e 2 and higher. 

Our starting point is the expression displayed in Sec. 5.1.5 [Eq. (413)] of LRR [4], 

= x')V a r - -J—U(x, x")V Q r adv + ± fo a U(x, x') + u a 'v a ,U(x, x')V a u 

2r z 2r adv 2r L 



V a U(x, x") + u a " V a »U{x, x")V a v\ + | [v(x, x')V a u - V(x, x")V a v 



-lj\ a V(x,z)dT. (2.1) 

We have introduced a large number of symbols. To begin, x is the field point at which the singular field is evaluated, 
and the world line is described by parametric relations z^{t) involving the proper-time parameter r. The points x' 
and x" on the world line are known respectively as the retarded and advanced points associated with x; these are 
defined such that x and x' are linked by a unique future-directed null geodesic originating on the world line, while x 
and x" are linked by a past-directed null geodesic that also originates on the world line. We define the scalar field 
u(x) as the value of the proper-time parameter at x' = z(t = u); this is known as the retarded time function of the 
field point x. Similarly, we define the advanced time function v(x) as the proper time at x" = z(t — ii). With o~(x, z) 
denoting Synge's world function [36], equal to half the squared geodesic distance between the field point x and the 
point z on the world line, we have that <j(x,x') = a(x,x") = 0. (Our subsequent developments rely heavily on a 
working knowledge of the general theory of bitensors; this material is reviewed in Sec. 2 of LRR [4].) 

The gradient of Synge's world function is denoted a a (x,x) if a(x,x) is differentiated with respect to its first 
argument, and cr a (x,x) if it is differentiated instead with respect to its second argument. The retarded distance r(x) 
between x and the world line refers to the retarded point x' and is defined by 

r := o- a ,{x,x')u a ' , (2.2) 

where u a is the particle's velocity vector at x'; this is an affine-parameter distance along the null geodesic that links 
x to x' . The advanced distance r a d v (x) between x and the world line refers instead to the advanced point x" and is 
defined by 

r-adv := -o- a "(x,x")u a , (2.3) 
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where u a is the particle's velocity vector at x"; this also is an affine-parameter distance. We note that the distance 
functions are both nonnegative. As a consequence of their defining relations (see Sees. 3.3.3 and 3.4.4 of LRR [4] for 
an extended discussion), we have the scaling relations r = 0(e), r a d v = 0(e), and v — u = 0(e). It also follows from 
the defining relations that the gradients of u, v, r, and r a d v that appear in Eq. (2.1) are given by 



V q m = -a a (x,x')/r, 

W a v = <r a (x,x")/r adv , 

V a r = a a ifjiu a V a u + o a i a u a , 

V Q r adv = -a a »/3»u a u fi V a v - a a » a '> 



(2.4) 
(2.5) 
(2.6) 
(2.7) 



Equations (2.4) and (2.5) are always valid and follow directly from the conditions a(x,x') = a(x,x") = 0; Eqs. (2.6) 
and (2.7) are valid when the world line is a geodesic of the curved spacetime. 

The biscalars U(x,z) and V(x,z) that appear in Eq. (2.1) are respectively the "direct" and "tail" parts of the 
retarded Green's function G(x,z) associated with the scalar potential $(x). (The general theory of scalar Green's 
functions in curved spacetime is reviewed in Sec. 4.3 of LRR [4].) For our purposes in this section, the only relevant 
properties of these objects are the scaling relations 

(2.8) 
(2.9) 
(2.10) 
(2.11) 
(2.12) 
(2.13) 
(2.14) 
(2.15) 
(2.16) 



U(x, x') 


= l + 0(e 4 


U(x,x") 


= l + 0(e 4 


V a U(x,x') 


= 0(e% 


V a 'U(x 1 x') 


= 0(e 3 ), 


V a U(x,x") 


= 0(e% 


\7 a „U(x,x") 


= 0(e% 


V(x,x') 


= 0(e% 


V(x,x") 


= 0(e% 


V a V(x,z) 


= 0(e). 



These relations are valid in a Ricci-flat spacetime only. 

Substituting Eqs. (2.8)-(2.16) into Eq. (2.1) produces the substantial simplification 

*a(*) = - ^Y V ^adv + 0(e 2 ). 



(2-17) 



This will be turned into a more explicit expression in the course of the following subsections. 



B. Reference point on the world line 

The expression of Eq. (2.17) refers to two separate points on the world line, the retarded point x' = z(u) and the 
advanced point x" = z(v); each is linked to x by the null conditions a(x,x') — a(x,x") = 0. We find it convenient to 
introduce a third point x = z(f) on the world line, and to go through the lengthy procedure of re-expressing V a $> s (x) 
solely in terms of tensorial quantities that are evaluated at x. An important aspect of this transcription is that we 
take the point x to be completely arbitrary, except for the following restriction: We assume that x is in a spacelike 
relation with x, so that it lies after x' but before x" on the world line; x is otherwise arbitrary. The transcription 
produces two major advantages: First, it consolidates the dependence of the singular field on the world line to a single 
point instead of two; and second, it eliminates the dependence of $^ on x that is only implicitly contained in x'(x) 
and x"(x). Our resulting expression for the singular field, which appears in Eq. (2.62) below, contains a dependence 
on x that is fully explicit. 

Having made arbitrary choices for the field point x and the reference point x on the world line, we define the 
quantities 

F :— (Tq(x, x)u & (2.18) 

and 

s 2 := (g & P + u & uP)o- a (x,x)a- p (x 1 x). (2.19) 
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We note that f = 0(e) and that its definition is similar to that of r and r a( jv provided in the preceding subsection. 
In fact, when x — > x' we have that f —* r, while f — > — r a( j v when x — ► x"; somewhere between x' and x" we have f 
changing sign, from a positive value to a negative value. The quantity s 2 = 0(e 2 ) is the squared distance between x 
and x as measured by an observer at x that is momentarily comoving with the charged particle. We note the useful 
identity s 2 — 2ct(x,x) + f 2 , which shows that s 2 is necessarily positive when x and x are in a spacelike relation. 

Our remaining task is to carry out the procedure described in the first paragraph of this subsection. This involves 
many steps, and lengthy calculations. The first step is to determine the positions of the retarded and advanced points 
relative to the reference point x; this we do in Sec. II C. The next steps involve computing the various pieces of the 
singular field of Eq. (2.17) and expressing them in terms of tensorial quantities defined at x. We begin in Sec. II D 
with a computation of r and r a( j v . We continue in Sec. II E with a computation of a a (x,x') and a a (x,x"), which 
appear in the expressions for V Q u and V a v — see Eqs. (2.4) and (2.5). In Sees. II F and G we compute the quantities 
a aa 'U a , a aa nu a , G a ipiu a , and a a "f3"U a vP , which are involved in the expressions for V Q r and V a^adv See 
Eqs. (2.6) and (2.7). In Sees. II H and I we collect our results and compute V Q u, V Q u, V Q r, and V a r a d v - And finally, 
in Sec. II J we produce our final expression for V Q $ s (x) — see Eq. (2.62) below. The reader who does not wish to 
go through these computations may simply jump to Sec. II J for the punch line. 



C. Retarded and advanced points 

We wish to determine the positions of the retarded point x' = z(u) and the advanced point x" = z(v) relative to 
the arbitrary reference point x = z(f) on the world line. We will achieve this by obtaining expressions for 

A + :=u-f>0 (2.20) 

and 

A_:=u-f<0. (2.21) 

These will be given in the form of expansions in powers of e. We will rely on Taylor-expansion techniques reviewed 
in Sec. 3.4 of LRR [4], as well as the standard bitensorial expansion (see, for example, Ref. [39]) 

^ - m - l R m^° 6 + ^mm-^* 5 ^ + °( £4 )> ( 2 - 22 ) 

in which the metric, as well as the Riemann tensor and its covariant derivative, is evaluated at the reference point x. 
We keep x fixed and introduce the function 

a(r):=a(x,z(r)) (2.23) 

of the proper-time parameter r on the world line. We note the special values cr(u) = <r(v) = 0, and that <t(t) is 
positive in the interval u < r < v. We express <t(t) as a Taylor expansion around the reference point r = f, and 
evaluate it at r = w, which stands collectively for either u or v. With A := w — f (and therefore equal to cither A + 
or A_), the result is 

= a + &A + iaA 2 + icr (3) A 3 + -^er^A 4 + • • • , (2.24) 

where a :— a(f) and all derivatives of c(r), which are indicated by overdots or a number within brackets, are evaluated 
at t = f . The computation of the derivatives is simplified by the fact that the motion is geodesic. From Eq. (2.23) 
we have & = a & u a , a = a a pu a u^ , = a Bi ^u a u^u^ 1 , and = a a 0^gU a u^u^u s . To evaluate the first derivative 
we simply involve Eq. (2.18) and get a = f. For the second derivative we involve Eq. (2.22) and obtain 

<j=-l- ^Ruaua + ^Ruaua\a + 0(e 4 ), (2.25) 

where we have introduced the notation R ua ua '■= Rap,f3v ua < jflu ^ a ' u an d Ru<jua\a '■— Rafil5v-.\ ua ' J ^ u ^ aV many 
variants of this notation will appear below. To evaluate the third derivative we begin by differentiating Eq. (2.22) to 
obtain an expansion for cr a ^, which we then contract with the velocity vector. The result is 

<J {3) =-\Ruau,\u + 0(e i ). (2.26) 
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We proceed similarly for the fourth derivative and obtain 

<r (4) = <3(e 2 ). (2.27) 

Gathering the results, Eq. (2.24) becomes 

= a + f A - 1 + ii? uaM(T - ^R uaual ^j A 2 - -^i? UCTMCT |„A 3 + 0(e 6 ). (2.28) 

This equation must now be solved for A. 

We assume that A can be expressed as an expansion in powers of e, in the form 

A = Ai + A 2 + A 3 + A 4 + C>(e 5 ) (2.29) 

with A„ = 0(e n ). Substituting Eq. (2.29) into Eq. (2.28) and equating each coefficient of e™ to zero returns a hierarchy 
of equations to be solved. The first equation determines Ai: 

A 2 - 2fAi - 2a = 0. (2.30) 

In view of Eq. (2.19) we have that 2cr = s 2 — f 2 , and the two solutions to Eq. (2.30) are A^ = f + s and A^ = f — s. 
The remaining equations produce A 2 = 0, 

D A 2 

Aa - w^r ( } 

. _ (Ali? M(Ttl(T | M — R U aua\<r) Aj ('2 32) 

4 _ 24(f- Ai) ' 1 ' ' 

The series expansion for A is now determined to the required degree of accuracy. 
Collecting our results, we conclude that A± is given by 

A± = (f ± s) T ^r^Ruau* T [(f ± s)R uau „\ u - R uaua \ a ] + 0(e 5 ). (2.33) 

This determines the positions of the retarded point x' — z(u) — z(f + A_) and the advanced point x" = z(v) = 
z(f + A + ) relative to the reference point x = z(f) on the world line. We note that in Eq. (2.33), the first term on the 
right-hand side is of order e, the second term is of order e 3 , and the third term (involving the square brackets) is of 
order e 4 . We recall the notation introduced below Eq. (2.25): In Eq. (2.33), 

Ruaua := flaofci^uV (2.34) 

and 

Ruaua\a ■= ; ,,X""-"" ; T"T A . (2.35) 

The notation is unambiguous and easily adaptable to other projections of the Riemann tensor; for example, R uiyua \ u := 



and 



R a ^ j? .^u a a fl u f3 a''u 1 also appears in Eq. (2.33). 



D. Calculation of r and r adv 

We now wish to express r and r a d v in terms of tensorial quantities that arc evaluated at x. Once more our strategy 
is to perform a Taylor expansion around r = f. From Eqs. (2.2) and (2.23) we obtain r = &(u), which may be 
expanded as 

r = a + aA_ + ^ 3 »A 2 + A 3 _ + • • • , (2.36) 

where A_ = u — f and where all derivatives of <t(t) are evaluated at r = f . Similarly, Eqs. (2.3) and (2.23) give 
^adv = -&(v) and 

radv = -&- ^A+ - ^ A 2 + - ^ (4) A 3 + + • • • , (2.37) 
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where A + = v — f. We recall that & = f, and that expressions for the higher derivatives of <t(t) were obtained in 
Eqs. (2.25)-(2.27). Furthermore, Eq. (2.33) gives A± as an expansion in powers of e. Making these substitutions into 
Eqs. (2.36) and (2.37) produces, after some simplification, 



f 2 - s 2 



r = s — 



6s 



24s 



(f + s)R 

U<7Ua\<7 

+ 0(e 5 ) 



(2.38) 



and 



-2 2 

6s 



r + s 
24s 



(r + s)(r - 2s) J R„ (TU(T |„ - (r - s)^^ + (9(e 5 ). 



(2.39) 



We note that in Eqs. (2.38) and (2.39), the first term on the right-hand side is of order e, the second term is of order 
e 3 , and the third term (involving the square brackets) is of order e 4 . Notice also that the difference between r and 
^adv is of order e 4 . 



E. Calculation of a a (x,x') and a a (x,x") 

We continue to keep x fixed and introduce the vector-valued function 

cr Q (r) := <r a (x, z(t)) 



(2.40) 



on the world line; the vectorial index refers to the fixed point x, and a a (r) is a set of four scalar functions of the 
argument r. In terms of this we have a a (x,x') — cr a (u) and a a (x,x") — <T a (v), and we wish to express these in 
terms of tensorial quantities evaluated at x. Letting w stand for either u or v, and re-introducing A := w — f, Taylor 
expansion gives 



o- a (w) =* a + a a A + ^a a A 2 + l^A 3 + ^^ 4) A 4 + 



(2.41) 



where a a — c Q (f), and where all derivatives are evaluated at r = f. We shall evaluate each term on the right-hand 
side of Eq. (2.41). We rely on results obtained in Sec. II C, as well as the standard bitensorial expansion (see, for 
example, Ref. [39]) 



(2.42) 



where g^Ax,x) is the parallel propagator, which takes a vector at x and carries it to x by parallel transport. 

We begin by recalling the identity a a (x, x) — —g a a (x, x)a a (x, x), which follows from the geometrical interpretation 
of a a and — cr a as tangent vectors on the spacelike geodesic that links the points x and x. The identity allows us to 
write a a as 



c Q = —g a cr & . 

The derivative of a a (r) is given by & a = (J aa u a , and involving Eq. (2.42) produces 



O-a = 9 , 



(2.43) 



(2.44) 



where we have introduced a notation similar to that of Eqs. (2.34) and (2.35): Rq 

R. 



R ap,p»< jflu ^°~ l/ and R&aua\o 



aftl3v;\ a u & & 

The second derivative of a a {r) is <j a — a aa pU a uP , and an expression for o~ a& p can be obtained by differentiating 
Eq. (2.22) with respect to x a . This expression is simplified with the help of Eq. (2.42), which can be truncated to 
a^ a = —g^a + 0(e 2 ) for the purposes of this computation. The end result, after contracting with u a u" and invoking 
the symmetries of the Riemann tensor, is 



(2.45) 
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Similar computations for the third and fourth derivatives produce 

1 



= q a 
"a ii a 



R(xucru\u 0(e 



and 



0(e). 



(2.46) 



(2.47) 



Equations (2.43)-(2.47) can now be incorporated into Eq. (2.41), in which we also substitute Eq. (2.33). After 
simplification, the final results are 



a a (x,x') 



<J<x + {r - s)u 6 



1 (f - s) 2 1 _ 2 

— (r S)Ra.cru(7 ~t~ T, Ruaua^a 77 S ) R*d 

6 6s 3 



±lf-8)R- , 
12 1 24s 



12 



(f-s) 3 i? e 



+ 0(e 5 ) 



(2.48) 



and 



a a (x,x") = g a a \ - <7 a + (r + s)u & 
+ 



1 (f + s) 2 1 2 

g if ~l~ s)Raau(T Ruaucr u a 7t(^ ^ ^) Rq 

| (f -|- s ) 2 / \ 

— (r + s)R &auiJ \ a H ^(r + s)i? M(Ttl(T | u — R uaua \ a J tia 



24s 



1 



1 



2^ ( r "I" s ) (3-RaMo-ti|cr ~t~ Raaua\u) + 1r) ( r + s ) Rauau\u 



12 



0(e 5 ) . 



(2.49) 



In these equations, terms grouped within square brackets are of the same order of magnitude: The first group of terms 
is of order e, the second group is of order e 3 , and the third group is of order e 4 . We recall the notation introduced below 
Eq. (2.44): In Eqs. (2.48) and (2.49), the various partial projections of the Riemann tensor are given by equations of 
the form 



and 



R&aua\a '■— R&fipv; U 13 (J V (J X . 



(2.50) 
(2.51) 



The full projections of the Riemann tensor have already been introduced in Eqs. (2.34) and (2.35). 



Calculation of a„„/u a and a„„nu a 



In terms of the vector-valued function a a (r) introduced in Eq. (2.40), we have that a aa 'U a = <J a {u) and a aa >>u° 
<j a (v). With w standing for either u or v, and with A := w — f , Taylor expansion gives 



& a (w) =a a + a a A + A 2 + \a<$ A 3 + • • • , 
2 o 



(2.52) 



in which all derivatives of a a (r) are evaluated at r = r. These quantities are displayed in Eqs. (2.44)-(2.47), and 
Eq. (2.33) provides an expression for A. Substitution and simplification yields 



— Rn 



^(r-s)R 6 



12 



a.aua\a 



-(f-a)(3ife 



■a.uau\<7 ~l~ R-acrur 



\u) ~t~ ^ (T s ) ^a.uau\u 



(2.53) 
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and 



a niy "U 



1 2 

~^ Raaua 77 (T s)Rq 
3 



12* 



Q:crU(7|(7 



+ 0(6 4 ) 



(2.54) 



In these equations, the first term on the right-hand side is of order e°, the second group of terms is of order e 2 , and 
the third group is of order e 3 . 



G. Calculation of <j„iaiu a u 13 and o„unnu a u 13 



Returning to the function ct(t) defined by Eq. (2.23), we have that a a > p>u a 'u 13 ' = <r{u) and a a » p»u a " u 13 " = a(v). 
Taylor expansion gives 



a(w) =a + a (3) A + ^ (4) A 2 + • • • . 



Substitution of Eqs. (2.33) and (2.25)-(2.27) produces 



<Ta'l3'U a U P = -1 - ^R U aua + Ruau<j\a ~ 3(r - s)R uau(T \ u + (3(e 4 ) 



and 



// of/ 
a a tffjiiu u p 



3 12 



- 3{r + s)R urTua \ u +0(e 4 ). 



(2.55) 



(2.56) 



(2.57) 



In these equations, the first term on the right-hand side is of order e°, the second term is of order e 2 , and the bracketed 
terms arc of order e 3 . 



H. Calculation of V Q w and V u 



According to Eq. (2.4), the gradient of the retarded time function u(x) is V a u = —<7 a (x,x')/r, and according to 
Eq. (2.5), the gradient of the advanced time function v(x) is \7 a v = a a (x, x")/r ac iy Expansions of r and r ac j v in 
powers of e were obtained in Sec. II D and presented in Eqs. (2.38) and (2.39); these can easily be converted into 
expansions for 1/r and l/r a d v - Expansions for <j a (x,x r ) and a a (x,x") were developed in Sec. II E and presented in 
Eqs. (2.48) and (2.49). Combining these expansions produces 



and 



s 9 



+ 



<r a + (r - s)u 
(f - s) 2 (f + 2s) 



+ 



g [v s) R a aua ^ (^* S ) ^auau "i" 



-2 2 

r A — s 
6s 2 



R 



6s 2 



RufTua^a 



Y2 ( r s )Rnou<j\cj + fi ( r s ) Rau<ju\o + 0/1 ( r s ) Raaua\u 



1 



1 

24 1 



(f-s) 



24s 2 



((f-s)(f + 2s).R. 



(f + s)i? u 



|<7 I "a 



+ r 24g 2 ((r - s)(r + 3s)R uaua]u - (r + 2s)i? UCTU<T | CT )uo 



+ 0(e 5 ) 



(2.58) 



\7 a v 



—9 



+ 



(Ta + (r + s)u. 

R 



+ 



g(r + s)i? 5 



{t -\- s) Rauau ~t" 



r 2 - s 2 



6s 2 



R 



(f + s) 2 (r -2s) 



6s 2 



i 



j y _|_ g / \ 

- j^(r + sfR &uau]u + ^^((f + s){f - 2s)R uaU(y \ u - (f - s)R uauala ja 6 



+ (if + s)(f- 3s)R U(TU „ lu - (f - 2s)R 



U(7ua\(T ) 



0(e 5 



(2.59) 
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Once more we have grouped terms of the same order of magnitude. The first group of (square-bracketed) terms within 
the curly brackets is of order e, the second group is of order e 3 , and the third group is of order e 4 . Noticing the factor 
,s _1 in front of the curly brackets, this means that V Q w and V a v contain terms of order e°, e 2 , and e 3 ; the neglected 
terms are 0(e 4 ). 



I. Calculation of V Q r and V a r adv 



The gradients of the retarded and advanced distance functions were defined in Eqs. (2.6) and (2.7); for example, 
V Q r = a a ifjiu a u 13 V a u+a a ' a u a . Each piece of this expression was calculated separately in the preceding subsections: 
<7 a 'l3'U a uP was computed in Sec. II G, V a w was computed in Sec. II H, and a a i a u a was computed in Sec. II F. 
Combining all these results, we arrive at 



V a r = 



-V 



+ 
+ 



12 



TfRn 



r 2 + s 2 



^TRaaucr ^{v S )R a ucru ~t~ ,. 



t{t 2 — s 2 ) 

Ruaua&a ~t~ g^2 Ruaua^c 



a + -(f 2 - S 2 )R, 



a.U(7u\a 



24 



24~s 2 ~ S ^ r2 +rs + As2 ) R u<?uv\u - (r 2 + s 2 )R uaua \ a ^a a 



+ ^| ((r ~ «)(^ 2 + 2fs + 3s 2 )i? M(Ttlff |„ - f(f + s)A 



■uaucr\a 



+ 0(e 5 ) 



(2.60) 



and 



-V. 



a a + ru a 



,7? -l(f 2 -s 2 )R- I ^ +s2 j? cr- + ^ lii? 

' J ^acri/o" V ^ J 1 l olu<ju \ n 2 ^ucrua u a \ n 9 J ^uaua u, c 



6s 2 



6s 2 



+ 



12 



-f7? 



1 



1 



cteru crj 



24 



j^O + s)(r 2 -rs + As 2 )R U(yua \ u - (r 2 + s 2 )i? M(Ttl(T | CT ) <J 5 



24s 
f + s 



+ t^j ((»" + s)(r 2 - 2rs + 3s 2 )R U(Jualu - r(r - s)R uauiJ \ a 



12 



+ 0(e°) . 



a.uau\u 



(2.61) 



In these equations, the first group of (square-bracketed) terms within the curly brackets is of order e, the second group 
is of order e 3 , and the third group is of order e 4 . Noticing the common factor s _1 , this means that V Q r and V a r a d v 
contain terms of order e°, e 2 , and e 3 ; the neglected terms are 0(e 4 ). 



J. Final result: The singular field 



Our final expression for the singular field is obtained by substituting the expansions of Eqs. (2.38), (2.39), (2.60), 
and (2.61) into Eq. (2.17). After a long computation and much simplification, we obtain 



KM 



2_ & 

B 3^ a 



o & + ru e 



+ 



1 1 2 2 3f 2 -S 2 

— T R(x aU(y ^ 3 jRaucru ~t~ 



r(r 2 — s 2 ) 

& + - 7. L R 



2s 2 



UGUG Lb Qt 



+ 



fR 

12 



QL<JUa\(J 



+ ^{r 2 - S 2 )R & uau\a + ^(^ 2 - S 2 )R &aU(T \ u - ~ 3s 2 )B\ 



auaulu 



-2 2 



+ 0(6 5 ) 



(2.62) 



In this equation, terms grouped within square brackets are of the same order of magnitude. The first group of terms 
is of order e, the second group is of order e 3 , and the third group is of order e 4 . Noticing the common factor s -3 , this 
means that $^ contains terms of order e -2 , e°, and e 1 ; the neglected terms are 0(e 2 ). 
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The dependence of <f>^ := V Q $ S on the field point x is contained in the common factor g a a (x,x), and also within 
the many occurrences of o~ & (x,x). This quantity appears explicitly in Eq. (2.62), and it is involved in the definitions 
of 

f := a & u a 

and 

s 2 := {grt + tfu^aaap, 

which were introduced in Eqs. (2.18) and (2.19), respectively. The gradient of the world function is also involved 
in the various projections of the Riemann tensor (and its covariant derivative) that appear in Eq. (2.62). We recall 
the notation introduced in the preceding subsections: A subscript u indicates a projection along u a , a subscript a 
indicates a projection along cr a (x, x), and a vertical bar indicates that the projection involves the covariant derivative 
of the Riemann tensor. For example, 

and 

The Riemann tensor, its derivatives, and the particle's velocity vector u a are all evaluated at the reference point x, on 
the world line. We recall that x is chosen to be in a spacelike relation with x, but that it is otherwise arbitrary. 

III. COORDINATE EXPANSIONS OF BITENSORS 

The calculations presented in Sec. II culminated into an expansion of the singular field 3>„(x) in powers of e, the 
distance between the field point x and the reference point x on the world line. This expansion is fully covariant, and 
the dependence on x is explicitly contained in cr a (x, x) and g a a (x, x). Our task in this section is to develop coordinate 
expansions for these two bitensors, in powers of 

w a :=x a -x a , (3.1) 

the difference in the coordinate positions of the points x and x. These expansions, of course, will not be covariant; they 
will depend on the choice of coordinate system. When the coordinate expansions to be obtained here are substituted 
into Eq. (2.62) for the singular field, the result will be an explicit expression for the expanded 3>q(x) m the adopted 
system of coordinates. 

A. Description of the geodesic linking x to x 

We assume that there exists a unique geodesic segment that begins at x, and ends at x. This geodesic segment is 
denoted /3, and it is described, in the adopted coordinate system, by the parametric relations p a (X). We assume that 
the parameter A is an affine parameter, and that it is limited to the interval < A < 1. We have that p a (0) = x a and 
p a (l) — x a , where x a are the coordinates assigned to x, while x a are the coordinates assigned to x. 

The functions p a (X) may be expressed as Taylor expansions about A = 0: 

p a {\) = p a (0) + P a (0)\ + \p a (0)\ 2 + ^ 3 >(0)A 3 + ^ a(4) (0)A 4 + • • • , (3.2) 

in which overdots, or a number within brackets, indicate repeated differentiation with respect to A. Equation (3.2) 
implies 

P a (X) =p a (0) +p a (0)X + V (3) (0)A 2 + V (4) (0)A 3 + • • • (3.3) 

and 



p«(A) = p a (0) + p a(3 \0)\ + -p«( 4 )(0)A 2 + • 



(3.4) 
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These quantities are linked by the geodesic equation, 

p a (A)+r« /37 (A)/(A)^(A) = 0, (3.5) 

in which T Q ^ 7 (A) is the Christoffel connection evaluated on (3. 

The range of the affine parameter A was chosen to ensure that the tangent vector p a (X) is intimately related to the 
gradient of Synge's world function. In fact, as reviewed in Sec. 2.1.3 of LRR [4] — see in particular Eqs. (55) and 
(56) — we have that a a (x,x) = g a p{x)p a (l) and 

cr a (x,x) = -g a p(x)p a (0). (3.6) 
B. Calculational strategy 

Our first goal in this section is to obtain an expansion of (Ja(x, x) in powers of the coordinate difference of Eq. (3.1). 
This is given by w a — p a (l) — p a (0), or 

w a = P a (0) + |p"(0) + ip tt(3) (0) + ^P Q(4) (0) + • • • • (3.7) 

We will achieve this in four steps. First (Sec. II C), we substitute Eqs. (3.3) and (3.4) into Eq. (3.5) and solve for 
p Q (0), p a ( 3 \0), and p a ^ 4 \0) in terms of p Q (0). Second (Sec. II D), we incorporate these results into Eq. (3.7) and 
obtain w a as an expansion in powers of p a (0). Third (also Sec. II D), we invert this series to obtain p a (0) expanded in 
powers of w a . Finally (Sec. II E), we substitute the result into Eq. (3.6); our final expression for a a (x, x) is displayed 
in Eq. (3.19) below. 

Our second goal in this section is to obtain an expression for g a a (x, x) expanded in powers of w a . This calculation is 
carried out in Sees. II F and II G and it follows a very similar strategy; our final expression for the parallel propagator 
is displayed in Eq. (3.30) below. 

The calculations that follow rely on the expansions displayed in Eqs. (3.2)-(3.4) and (3.7). We also will need the 
Taylor expansion of T^ 7 (A) about A = 0. Starting with 

r^ 7 (A) = r^ 7 (o) + ry^(o) [^(A) - ^(o)] + V^o) [ P »(\) - p"(o)] [^(A) - P "(o)] + • • • 

and involving Eq. (3.2), we arrive at 

r^ 7 (A) = r^o) + [1^(0)^(0)] a + \ [^^(0)^(0)^(0) + r a 0J jo)r(o)} a 2 + • • • , (3.8) 

in which the connection and its derivatives are evaluated at x on the right-hand side of the equation. 

C. Calculation of p a (0), p a(3> (0), and p a(4) (0) 

We substitute Eqs. (3.3), (3.4), and (3.8) into Eq. (3.5) and collect terms that share the same power of A. Setting 
the coefficient of the A term to zero yields the condition = p a (0) + r Q ^ 7 (0)p /3 (0)p 7 (0), or 

p Q (0) = -r^ 7 /(0)p 7 (0). (3.9) 

It is understood that here, the connection is evaluated at the reference point x. The same comment will apply below 
to all derivatives of the connection. 

Setting the coefficient of the A 1 term to zero yields 

- p^ 3 \0) + 2T a y (0)iF(0) + (0)p"(0)p"(0), 

and involving Eq. (3.9) gives 

p a(3) (0) = _ r ^ 75 j^(0)pT(0)p«(0) (3.10) 

with 

r^^r^^rvV (3.n) 
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Setting the coefficient of the A 2 term to zero yields 

o = p«( 4 )(o) + 2r" /37 [/(o)^( 3 )(o) + /(o)^(o)] +r /Jw [4/(o)^(oK(o)+/(o)f(o)f(o)] 

+ 1-^^(0)^(0)^(0)^(0), 
and involving Eqs. (3.9)-(3.11) gives 

p a(4) (0) = _r^ e /(0)^(0)p 5 (0)p e (0) (3.12) 

with 

pa pa /ipa FA 4 F Q FA 4 OF Q FA 1 _L AV a FA 1 F^ _L OF Q FA 1 F" 1 Q"\ 

1 /3 7 <5e •— 1 /3 7 ,5e — 41 /3m,7 <5e — 1 07, M Se ~ Zi 7 5,c + 41 0A* 7^ <5e + Zi [iv 0y L <5c 

It should be noted that ^ a p 1 s and 5e , as defined by Eqs. (3.10) and (3.12), are both fully symmetric in their 
lower indices. This symmetry has not, however, been implemented on the right-hand sides of Eqs. (3.11) and (3.13). 
Although this could easily be achieved, this operation is not necessary and we opt to leave these expressions as they 
are. 



D. Calculation of p a (0) 
Combining Eqs. (3.7), (3.9), (3.10), and (3.12) gives 

™ Q = p a (o) - \r%f(o)p~<(o) lry^(o)^(o)/(o) - ^O 75e /(o)^(o)/(o)p e (o) + • • • . (3.w) 

This is an expansion of w a in powers of p a (0). The inverted series will take the form of 

p°(0) = w a + A^wPvf + A^w^w^w 5 + A a pi&t w w^w s w e + • • • , (3.15) 

and the coefficients A a p , A a p lS , and A a p lSf can be determined by inserting Eq. (3.15) into Eq. (3.14) and demanding 
that the substitution returns the identity p a (0) = p a (0) . 

Elimination of the quadratic terms gives rise to the condition 

A a 01 := \t% t (3.16) 

Elimination of the cubic terms produces A a p lS = \T a p lS + A a p li T ,i s . This becomes 

after involving Eqs. (3.11) and (3.16). Elimination of the quartic terms produces 

Act 1 pa , 1 a a pA 1 ^ Aa -pH ni/ . ^ \ct pA 1 i ^ Act pA 1 i ^ Act pA 1 

A f3y8e ~ 24 1 /3 7 <5e + 3 0f* 7<5e ^ /3-y 1 5e + 2 |87|U 5e 2 ^ 7<5 2 ^ 5e ' 

This becomes 

^ Q /3 7 «5 e : = ( rQ /37,<5e + Ol,/^ + ^"flA^V.e + ^"n^^-y^Se) (3.18) 

after involving Eqs. (3.11), (3.13), (3.16), and (3.17). 

Equation (3.15), with the coefficients of Eqs. (3.16)-(3.18), gives the expansion of p a (0) in powers of w a ~ x a — x a . 
We recall that the coefficients involve the Christoffel connection and its partial derivatives evaluated at the reference 
point x. We also remark that while A a p s and A a p Se have been defined in Eq. (3.15) as being fully symmetric in 
their lower indices, the right-hand sides of Eqs. (3.17) and (3.18) have been left in a non-symmetric form. 



E. Final result: cr^ expanded in powers of w a 

The gradient of Synge's function is obtained by substituting Eq. (3.15) into Eq. (3.6). The result is 

— cr g (x, x) — gapvjP + AapyW^w 7 + Aa^sw^w 1 w 5 + A a pys e w^ l w n 'u/ ' w e H , (3.19) 

where the coefficients A a p J} A a pjs, and A a pjSe are obtained from Eqs. (3.16)-(3.18) by lowering the first index with 
the spacetime metric (as if these quantities were tensors). We recall that the metric g a p, as well as the connection 
r Q ^ 7 and its partial derivatives, is evaluated at the reference point x. Equation (3.19) is the required coordinate 
expansion of <7q in powers of w a — x a — x a , the difference in the coordinate positions of the points x and x. 
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F. Parallel transport on the spacelike geodesic 

We next turn to our second task, the development of a coordinate expansion for the parallel propagator g a a (x,x). 
We begin by introducing an arbitrary dual vector q a (X) that we take to parallel transported on (3, the geodesic segment 
that links x to x. The definition of the parallel propagator implies that q a (l) = q a (x) and q a {0) = q&{x) are related 
by 

q a (x) = g a a (x,x)q a (x). (3.20) 

We shall calculate the parallel propagator by expanding q a (X) in a Taylor series about A = 0, evaluating this at A = f , 
and comparing the result with Eq. (3.20). 
The Taylor expansion is 

9a(A) = g a (0) + q a (0)X + ^'"(0)A 2 + i#(0)A 3 + • • • (3.21) 

and it implies 

q a (X) = q a (0) + q a (0)\ + ^ 3) (0)A 2 + • • • . (3.22) 

The dual vector is parallel transported on f3 if 

q a {\) - r%(A)^(A)/(A) = 0. (3.23) 

We shall work on this equation, using the expansions for pP{\) and rAt Q/3 (A) that are displayed in Eqs. (3.3) and (3.8), 
respectively. 

We make the substitutions in Eq. (3.23) and collect terms that share the same power of A. Setting the coefficient 
of the A term to zero yields the condition 

fc(0)=r' t a/3 ^(0)/(0). (3.24) 

We recall that the Christoffel symbols are evaluated at the reference point x; the same remark applies to their 
derivatives, which will appear in expressions below. 
Setting the coefficient of the A 1 term to zero yields 

fe(o) = 1% Mo)/(o) + ^(o)/(o)] + rvx(o)^(o)/(o). 

This becomes 

9a(0) = g" a/j7 g M (0)^(0)^(0) (3.25) 

with 

Q^ap-f = _ ^^au^fSj + ^^-fv^afii (3.26) 

after involving Eqs. (3.9) and (3.24). 

Setting the coefficient of the A 2 term to zero yields 

li 3) (0) = r^[g M (0y< 3 >(0) + 2^0^ 

+ [rV/(°) + ^ ,^(0)^(0)] ^(o)/(o). 

This becomes 

<zi 3) (0) = Q% 7tf ^(0)^(0)^(0)^(0) (3.27) 

with 

QA 1 pM _ pM V v _i_ F^ 1 — OF" F^ _1_ orf pf ~pv pA 4 

a/3 7< 5 ~ 1 a/3, jS 1 cw 1 /3 7 ,<5 * 1 a/3 1 „ 7 ,a zi 0~f L av.S zl /31/ 1 a 7 ,5 1 /3-y 1 

I 9PM pz> pA npP p" pA pM pA , pM pf pA /q OQ\ 

+ lL av L (j\ L ~f6 ~ ZL f3v L a\ L 7 5 _ 1 U \ L afj L 7 5 + 1 /3,, 1 7 A i a<5> 
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after involving Eqs. (3.9)-(3.11), as well as Eqs. (3.24)-(3.26). 
Equations (3.21), (3.24), (3.25), and (3.27) combine to give 

?a(i) = g«(o) + r% ^(o)/(o) + ^(o)/(o)^(o) + \q^ 0i5 <^(o)/(o)^(o)/(o) + • • • . 

There is a common factor of (^(0), and this equation has the same form as Eq. (3.20). The parallel propagator is 
therefore identified as 

gP a {x,x) = n + T%/(0) + iQ% 7 /(0)^(0) + \cT a ^(0)f(0)p s (0) + ■ ■ ■ • (3.29) 
The coefficients of this expansion in powers of p a (0) are given by Eqs. (3.26) and (3.28). 

G. Final result: g^ a expanded in powers of w a 

Our final expression for the parallel propagator is obtained by substituting Eq. (3.15) into Eq. (3.29) and expanding 
the result in powers of w a . After a lengthy calculation that involves Eqs. (3.16), (3.17), and (3.18), we obtain 



g* a (x,x) = n + + B%/«,' + B» afhsW Pw< w 5 + ■■■ (3.30) 



with 



B% := (3.31) 
£% 7 == + (3-32) 

nC ^ f QpM I OV u F^ 1 V u F^ 1 _l_ AT^ T v _l_ F' y F^ 1 

D afiyS ■— 22 \ «/3,7<5 + a/3 ^7,<5 07 av,5 ' * L flv L aj.S ' 1 Pf L aS.v 

■pH -pv pA | -pfi -pv pA | r)p/J p^ pA \ (n no\ 

_i /Si/ 1 aA 1 jS + 1 i/A 1 a/3 1 7 5 + zi /3,, 1 7A 1 - a<5 J ■ (6.66) 

The Christoffel connection and its partial derivatives are all evaluated at the reference point x. 

IV. SPHERICAL-HARMONIC DECOMPOSITION OF THE SCALAR FIELD 

Our task in this section is to identify a useful way of relating a spherical-harmonic decomposition of the field 
$ Q := V a $ to a spherical-harmonic decomposition of the potential As we shall see, the selected relation involves 
a decomposition of the vector <I> Q in terms of a tetrad of orthonormal vectors e?-. . So instead of decomposing $ Q 
in a set of vectorial harmonics, we shall decompose each frame component := $ e? > of the vector field — a 
scalar function of the spacetimc coordinates — in scalar spherical harmonics. The decomposition of $ Q in vectorial 
harmonics would make a viable alternative strategy, but one which would prove less convenient for the purposes of 
calculating regularization parameters — see Sec. V and the Appendix. Our scheme leaves open the choice of tetrad, 
which is a priori arbitrary; our particular choice is motivated by a desire to keep the spherical-harmonic decompositions 
of $ and <&( M ) as closely linked as possible. 

A. Spherical- harmonic decompositions 

Let &(x a ,9 A ) be a scalar field on a spherically-symmetric spacetimc. The spacetime manifold has the product 
structure M 2 x S 2 , in which M 2 is a two-dimensional submanifold that is orthogonal to the two-spheres S 2 . We let 
x a stand for any coordinate system that charts an open domain of the submanifold M 2 ; the lower-case Latin index a 
runs from to 1. We let 9 A be angular coordinates on the two-spheres; the upper-case Latin index A runs from 2 to 
3. For a Schwarzschild spacetime charted with the usual coordinates [t, r, 9, <p\, we have x a = [t, r] and 9 A = [9, </>]. We 
shall leave the coordinates x a arbitrary for the time being, but we adopt the canonical angular coordinates 9 A = [6, <f>]. 

We suppose that the scalar field is expressed as a decomposition in spherical-harmonic functions Y lm , so that 

$>{x a , 9 A ) = & m {x a )Y lm {6 A ). (4.1) 
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The sum over the integer I extends from I = to I = oo, while the sum over the integer m ranges from m = —I 
to m = I. To keep $ real the mode functions $' m must satisfy $ l <~ m = (— l) m $' m ; i n which an ovcrbar indicates 
complex conjugation. The gradient V a <& of the scalar field possesses the components 

9„$ = d a <$> lm Y lm (4.2) 

Ira 

and 

d A <I> = J2^ lm d A Y lm . (4.3) 

Im 

These relations show that a natural basis of expansion for V Q <I> would involve the scalar harmonics Y lm in the 
M 2 sector, and the vectorial harmonics d A Y lm in the S 2 sector. Following this route, however, would introduce 
complications at a later stage (refer to the last paragraph of Sec. 2 in the Appendix), and we shall adopt an alternative 
strategy. 

We introduce a tetrad of orthonormal vectors e? ^ at every point in the spherically-symmetric spacetimc. The 
superscript a is the usual vectorial index, and the subscript (/x) = {(0), (1), (2), (3)} is a label that designates an 
individual member of the tetrad. These vectors satisfy 

9ape^ v) = r) M (v) (4.4) 

with T)(ft)(v) = diag[— 1, 1, 1, 1]. The four quantities 

are the frame components of the vector V a ^>; these are scalar functions of the spacetime coordinates. The vector field 
can be reconstructed from its frame components by involving the dual tetrad e^ , which is defined by 

e^:=r,^g a ^ v) , (4.6) 
where fjMM = diag[— 1, 1, 1, 1] is the matrix inverse of f?^)^)- It is easy to show that 

V Q * = (4.7) 

Because each frame component is a scalar quantity, it is natural to decompose it in scalar harmonics. We 
therefore write 

* w {aT,e A ) =Y t ${™ ) (x a )Y lm (6 A ), (4.8) 

Ira 

and seek to determine the relation between $[™), the modes of the frame components, and $' m , the modes of the 
original scalar field. 

The answer is provided by substituting Eq. (4.5) into the equations = J $( M )T' m dfl, where d£l — sm0d9d(j> is 
an element of solid angle. After involving Eqs. (4.2) and (4.3), we obtain 



I'm' 

where the coupling coefficients are given by 



*(S = E [Cfo(l'm'\lm)d a &' m ' + C M (l'm'\lm)^' m '] , (4.9) 



C^(l'm'\lm) := J e^Y 1 '™ Y lm dU (4.10) 



and 

rl'mf \rlm 



C M (l'm'\lm) := J ef^d A Y l m Y lm dSl. (4.11) 

Here, denotes the (0, 1) components of each basis vector, while represents the angular components. The 
computation of the coupling coefficients requires the specification of the tetrad. 

In the rest of this section we will make a specific choice of tetrad (Sec. II B), compute the coupling coefficients for 
this tetrad (Sees. II C and II D), and give an explicit form to Eq. (4.9). The tetrad is displayed in Eqs. (4.21)-(4.25) 
below, and the resulting relation between <j>|™^ and $ /m was already displayed in Eqs. (1.23)-(1.26). 
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B. Choice of tetrad 



The choice of tetrad is in principle free, but we wish to find a tetrad that leads to a simple structure for the coupling 
coefficients. Specializing to Schwarzschild spacetime and the usual coordinates [t,r,6,</>], a possible choice of tetrad 
would be the usual orthonormal frame 



[r 1/2 , 0,0,0], 

[o,/ 1/2 ,o,o], 

[0,0,r-\0], 
[0,0,0, (rsintf)- 1 ], 



where 



/=-i-2£ 

r 



(4.12) 
(4.13) 
(4.14) 
(4.15) 

(4.16) 



It is easy to show, however, that while this tetrad would lead to simple (diagonal) coupling coefficients Cfy, it would 
also lead to coefficients C( M ) that couple each (Im) mode to an infinite number of (I'm') modes. We shall not, therefore, 
make this choice of tetrad. 

We shall instead introduce a "Cartesian frame" that is linked to the "spherical frame" of Eqs. (4.12)-(4.15) by the 
same relations that would hold in flat spacetime. Explicitly, our choice of tetrad is 



'(0) 

'(1) 

'(2) 
'(3) 



sin 9 cos cj> e 
sin 9 sin 6 e 



(r) 
a 

(r) - 



COS V COS ( 



sin ( 



cos 9 sin <j> e" S ) + cos <j> e 



'(0)' 



cos W e 



(r) 



sin We 



(<?)• 



(4.17) 
(4.18) 
(4.19) 
(4.20) 



We may loosely think of e?U as pointing in the "a; direction," of e? 2 \ as pointing in the "y direction," and of e" 3 ^ 
as pointing in the u z direction," with [x,y, z] representing a quasi-Cartesian frame related in the usual way to the 
quasi-spherical coordinates [r, 9, (f\ . Independently of this heuristic interpretation, we note that the transformation of 
Eqs. (4.17)-(4.20) defines a legitimate set of orthonormal vectors. And as we shall see, this tetrad has the desirable 
property of leading to a simple structure for the coupling coefficients. 

Combining Eqs. (4.17)-(4.20) with Eqs. (4.12)^(4.15) produces the following explicit expressions for the basis 
vectors: 



(4.21) 
(4.22) 
(4.23) 
(4.24) 



a _ 

(o) - 


^f, 0,0,0 






a 

(1) - 


0,77 sin0 


cos 4 


>, - cos 6 
r 


a 

(2) - 


0, y/f sin 9 


sin^ 


, - cos 9 

r 


a 

(3) - 


0, v 7 / cos e 


1 

r 


sin 9, 



sin ( 



e (0 \ = 



r sin 

COS( 



rsinf 



It is useful to introduce, as substitutes for e"^ and e" 2 ^ , the complex combinations 

±ie ±i4 



„CK _|_ • Oi 

e (l) ± le (2) 



0, VTsintfe^-cosfle^ 



r sinf 



(4.25) 



In terms of the complex vectors we have = [e" + ^ + e"_- ) ]/2 = Re[e™ + j] and e™ 2 ) = [e" + ^ — e"_^]/(2i) = Im[e" + ^]. In 
the sequel we will work primarily in terms of the complex tetrad e" ^, e" + j, and e%y 
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C. Calculation of C a (ll) (l'm'\lm) 

We may now substitute the tetrad of Eqs. (4.21)-(4.25) into Eq. (4.10) and calculate C"^(l'm'\lm), the first set of 
coupling coefficients. Our computations will rely on the standard identities (see, for example, Sec. 12.9 of Ref. [40]) 



^ nft J4> v lm _ _ l (l + m+l)(l + m + 2) l+1 , m+1 / (Z - m)(l - m - 1) ,_ 1%m+1 , 
&my ° V (« + l)(2/ + 3) + V (2Z-1)(2Z + 1) Y ' (4 - 2 ° 



rin flp-^V 1 " 1 / (^-m+l)(i-m + 2) ; + i, m _i (l + m)(l + m- 1) ;_i, TO _i , . 

" V (2Z + l)(2Z + 3) y V (2Z-1)(2Z + 1) ' ( J 



involving spherical-harmonic functions. 

Substituting Eq. (4.21) into Eq. (4.10) and invoking the orthonormality relations of the spherical harmonics reveals 



that the only nonvanishing component of C? >.(l'm'\lm) is 



'(0) ^ 

Substituting Eq. (4.25) into Eq. (4.10) and involving Eq. (4.27) shows that 



Cor^fcU. (4.29) 



_ (l + m-l)(l + m) r- / (/-m+l)(/-m + 2) R 

- "V (2Z-1)(2Z + 1) V/*i',i-i^, m -i + (2Z + l)(2Z + 3) V/<W<W, m -i (4-30) 



(2Z-1)(2Z + 1) VJ ' ' y (2Z + 1)(2Z + 3) 

is the only nonvanishing component of C^(l'm'\lm). Similarly, 



CT_ 



/ (Z-ro-l)(i-m) ns / (Z + m+l)(Z + m + 2) ^ 

(") = V (2Z-1)(2Z + 1) VfSv.i-iSm'^i -y (2Z + l)(2Z + 3) V/^',m^, m+ i (4.31) 

is the only nonvanishing component of (7"_^(Z'm'|/m). Finally, substituting Eq. (4.24) into Eq. (4.10) and involving 
Eq. (4.26) reveals that 



/ (l-m)(l + m) r- (l-m+l)(l + m + l) r- 

°^ ~ V (2l-l)(2l + l) Vf5l '< l - lSm ' m + V (2Z + l)(2Z + 3) V/ft',i+i^ OT (4-32) 



is the only nonvanishing component of C^(Z'm'|Zm). 



D. Calculation of C( M ) (i'm'|/m) 

We now substitute the tetrad of Eqs. (4.21)-(4.25) into Eq. (4.11) and calculate C^(l'm'\lm), the second set of 
coupling coefficients. These computations also will rely on the standard identities listed in Eqs. (4.26)-(4.28), as well 
as (see, for example, Sec. 12.7 of Ref. [40]) 

e^ide + i cot 6d^)Y lm = y/(l - m)(l + m + l)Y l ' m+1 (4.33) 

and 

e _i *(0 fl - \coted <p )Y lm = -yJ{l + m){l-m+l)Y^ m -\ (4.34) 
Substituting Eq. (4.21) into Eq. (4.11) immediately gives 

C (0) = 0. (4.35) 
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Substituting Eq. (4.25) into Eq. (4.11) produces an integral for C( + ) that involves the factor 

cosOe^deY 1 m + -^dsY 1 m . 

sva.6 

The first term can be expressed as 

e^de (cos 9Y l ' m ') + B mOe i +Y l ' m ' . 
The second term, on the other hand, can be expressed as 

ie i<#> cot 6»a (cos 6Y l ' m ') + i sin Bc^d^Y 1 '" 1 ' . 

The sum becomes 

e^idg + looted^) (cos BY 1 ' m ') - (to' - 1) sin 6»c i0 y'' m ' , 

and this is in such a form that Eqs. (4.26), (4.27), and (4.33) can now be involved. Multiplying by Y lm and evaluating 
the integrals returns 



(Z + m -l)(Z + m ) Z-l (/_ m+ i)(;- m + 2 ) Z + 2 
= V (2Z-1)(2Z + 1) — ^.'-^-'."-i + V (2Z + l)(2Z + 3) ~ ^ ( 4 ' 36 ) 



We similarly obtain 



" "V (2Z-1)(2/ + 1) ~ -y (2Z + l)(2Z + 3) — ( 4 - 37 ) 

Substituting Eq. (4.24) into Eq. (4.11) produces an integral that involves the factor 

sm6d e Y l ' m ' = dg(sm9Y l ' m ') - cos 9Y V w . 

This can be expressed as 

e _i *(aj - icot6»a )(sin6»e i,#, r z ' m ') +ie-^cot^(sin0e^y'' m ') - cos6»F'' m ', 

or as 

^(de - i cot (sin fle'^y'' m ') - (m' + 2) cos tfy''" 1 ' , 

which is now in a useful form. After involving Eqs. (4.26), (4.27), and (4.34), then multiplying by Y lm , and finally 
evaluating the integrals, we arrive at 



(l-m)(l + m) l-l (l-m+m + m+l) l + 2 
C(3) -~V (2«-l)(2/ + l) r ^'.'-i^ + y (2Z + l)(2/ + 3) ~ <W<^'- ( 4 - 38 ) 



E. Final result: in terms of $' m 

We substitute Eqs. (4.29)-(4.32), (4.35)-(4.38) into Eq. (4.9), which gives <j>|™ )7 the spherical-harmonic modes of 

the frame components $( At ), in terms of <I>' m , the modes of the original scalar field After evaluating the sums over 
I' and to', we find that the relationship is given by Eqs. (1.23)— (1.26), which are displayed back in Sec. I E. Inspection 
of these equations reveals that the relationship is simple: The structure of the coupling coefficients is such that <&'™j 

is linked to the neighboring modes $' ±1,m and <fr l±1 > m±1 only. This simplicity is a benefit of the choice of tetrad made 
in Sec. IV B; as was discussed in that subsection, other choices would produce more complicated relationships. 
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V. REGULARIZATION PARAMETERS 



The self-force acting on a scalar charge q moving on a geodesic of the Schwarzschild spacetime is proportional to 
'■— V a $ R , where <j> R is the regular potential that remains after the singular potential $ s is subtracted from the 
retarded potential <J>. A local covariant expansion for $^ := V Q $ S was worked out in Sec. II, and with the help of the 
results presented in Sec. Ill, this can be turned into an explicit expansion in Schwarzschild coordinates. In Sec. IV 
we introduced a tetrad of orthonormal vectors e?-, to resolve the vector fields $ Q , <I>|, and <I> R in terms of their frame 

components <&(^) := $ a e™^, <I>^ := 3> R e^, and := $|e"^, respectively. We have 

Also in Sec. IV we showed how the spherical-harmonic modes of the frame components can be related to those of 
the scalar potential Our task in this section is to compute the spherical-harmonic modes of the singular field 
and to extract from them the quantities known as regularization parameters. We will rely on the results obtained in 
Sec. II, III, and IV, as well as multipole-decomposition techniques reviewed in the Appendix. Most of the computations 
that are described below were carried out with the symbolic manipulation software Maple, with the help of the tensor 
package GRTensorII [37]. We will describe how these calculations were performed, but space considerations compel 
us to leave most details hidden. 



A. Definition of the regularization parameters 

The scalar charge q moves on an arbitrary geodesic of the Schwarzschild spacetime. We place this geodesic in the 
equatorial plane, and we assign to the particle the coordinates t = t(Y), r = r(r), 6 = ^ , and <f> = ip(r), in which r is 
proper time on the geodesic. These functions are determined by integrating the geodesic equations, which take the 
form t = E/f, r 2 = E 2 — f(l + L 2 /r 2 ), and ip = L/r 2 , in which an overdot indicates differentiation with respect to 
t. The constant E is the particle's conserved energy per unit rest-mass, and the constant L is the conserved angular 
momentum per unit rest-mass. We also have introduced the metric function / := 1 — 2M/r and its value f at r — r(r). 
At some instant r = To, the particle is found at the point x = [to, ro, f , 4>o] on its world line. We let /o := 1 — 2M/ro, 
and we wish to evaluate the scalar self- force at that instant. 

We are interested in the value of <I> R n at the point x. We calculate this by first decomposing the retarded field 
into spherical- harmonic modes, then subtracting the modes associated with the singular field and finally 

summing over all modes. This mode-sum converges because the regular field is smooth on the world line; the mode- 
sums associated with the retarded and singular fields do not converge on their own, because both fields diverge on the 
world line. In spite of this singular behavior, each spherical-harmonic mode $( M )/ m and $L)j m is bounded as X ^ X\ 
the singularity merely gives rise to a jump discontinuity of each mode at x. The value of each mode is therefore 
ambiguous at x, but this ambiguity is of no consequence because the discontinuity disappears when the mode of 
the singular field is subtracted from the mode of the retarded field. The modes of the regular field are continuous 
(and differentiable) at x because is smooth on the world line. In practice the discontinuity can be handled by 
evaluating each mode of the retarded and singular fields slightly away from x, performing the subtraction, and then 
taking the limit x — > x. 

A practical implementation of this prescription is contained in 

(5.2) 
where 

l 

®(M)r-= E *(^o,ro + A)y im (f,0 o ) (5.3) 

m— — I 

are the multi-pole coefficients of the retarded field, while 

*U ■= E *Wm(*o.ro + A)r Im (§,^) (5.4) 

m— — l 

are the multipolc coefficients of the singular field. The quantities $( A1 ); Tn and ^?^)j m are the spherical-harmonic modes 
of the retarded and singular fields, respectively. (Details regarding the spherical-harmonic decomposition of a scalar 



*f M) (*o, ro, f , <Ao) = Hm ]T fcf^ 
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function on S 2 are provided in the Appendix.) In Eqs. (5.2)-(5.4) we choose the displaced point to be at the same 
time coordinate t as x, but at a displaced radius r + A and a displaced azimuthal angle 



cA; 



(5.5) 



the constant c will be selected for convenience below, in Sec. V D — see Eq. (5.20). (The idea of introducing a 
displacement in the cj> direction goes back to Mino, Nakano, and Sasaki [19].) 

Equations (5.2)-(5.4) are at the core of the procedure to calculate the self-force. We imagine that the modes < I'( AI ); m 
of the retarded field can be computed with a convenient numerical method. From these we obtain the multipolc 



coefficients $ 



from which we subtract ^r^p 



the multipole coefficients of the singular field. These can be computed 



analytically, which is our task in this section. As wc shall see, they have the form 



(I + \)A W + B 



a 



(m) 



00 



D 



(z + D 



(i-m+i) 



(5.6) 



in which the coefficients A 



00' B (p)> 



(p), and £) (M ), 



known as regularization parameters, are independent of I but 



depend on the state of motion of the particle at x. The subtraction produces ^f^i, which we sum over all values of I 



to get $P\ and eventually the self-force. 



B. Rotation of the angular coordinates 

Efficient techniques to compute multipole coefficients were devised by Detweiler, Messaritaki, and Whiting [22]; 
these are reviewed in the Appendix. They rely on a rotation of the angular coordinates that maps the special point 
(6 = t£, <f> = (f>' ) to the North pole of the new angular coordinates. (This idea goes back to Barack and Ori [17, 18].) 
This rotation is described by 

sin# cos((/> — 4> ) = cos a, (5-7) 
sin 8 sin(0 — <j)' Q ) = sin a cos 0, (5.8) 
cos 6 = sin a sin 0, (5.9) 

which is a slightly modified version of the rotation described by Eqs. (A.5)-(A.7). The new angles are a and 0, and 
it is easy to see that the rotation does indeed map the point (6 — 4>b) to (a = 0, =?), with undetermined. 

As reviewed in the Appendix, each multipole coefficient ^f^i is computed by first expressing the singular field 
as a function of the angles a and 0, then extracting the Legcndre projection 

i r 1 

■= 2^ + 1 ) J ^l ) (t^r o + A,a,0)P l (cosa)dcoaa, (5.10) 



and finally averaging over the angles 0, 



1 f 27T 

*UW)d0. (5.11) 



The decomposition of (to, tq + A, a, 0) in Legendre polynomials relies on the techniques reviewed in Sec. 4 of the 
Appendix. The averaging over all angles relies on the techniques reviewed in Sec. 5 of the Appendix. A summary 
of the key results is provided in Sec. 6 of the Appendix. 



C. Calculation of "^(to^o + A, a,/?) 

The starting point of these computations is the calculation of the singular field at a position x, to which wc 
assign the (unrotated) coordinates [to, r + A, 0, cf\. This point is slightly displaced from x, to which we have assigned 
the coordinates \po,ro, ^, 4>o]- The displacement vector is 

w a :=x a -x a = [0,A,0- f, <£-<&)], (5.12) 

and an expression for <fr? Jx) expanded in powers of w a can be obtained by combining the results displayed in Sec. II 
J [Eq. (2.62) and following], Sec. Ill E [Eq. (3.19)], Sec. Ill G [Eq. (3.30)], and Sec. IV B [Eqs. (4.21)-(4.25)]. The 
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computation of the singular field is extremely tedious, and was handled by the symbolic manipulator GRTensorII 
[37] operating under Maple. 

The singular field is now expressed in terms of w e and (in addition to A), but these components of the 
displacement vector are functions of a and (3 that can be determined from Eqs. (5.7)-(5.9). We have 

w e = — arcsin(sinoi sin/3) (5.13) 

and 



;4> 



sin a cos (3 
V \J\ — sin 2 a sin 2 (3 



W r = arcsm : — cA, (5-14) 



where the (as yet unassigned) constant c was introduced back in Eq. (5.5). These quantities are small wherever a 
and A arc small, and the expansion of in powers of w a is valid in a small neighborhood of the North pole 
(a = 0,(3=1). 
Defining 



Q := Vl — cos a, (5.15) 
we observe that w and can each be expressed as an expansion in powers of Q: 

w e = -V2Q sin/3 - ^jj<2 3 (l -4cos 2 /3) sin/3 + 0(Q 5 ), (5.16) 

/2 

= -cA + V2Qcos/3+^Q 3 (9-8cos 2 /3)cos/3 + 0(Q 5 ). (5.17) 

We note that Q is formally of the same order of magnitude as a. Making these substitutions within the singular 
field returns a double expansion in powers of A and Q, which are formally considered to be of the same order of 
magnitude. This new representation of involves the rotated angles a and (3, which is required for its substitution 

into Eqs. (5.10) and (5.11). Moreover, the singular field is expressed in terms of sin/3, cos/3, and Q, functions that 
are globally well-defined on the sphere. This property is critical for the successful decomposition of in Legendre 
polynomials, and the subsequent averaging over (3. We stole this powerful idea from Detweiler, Messaritaki, and 
Whiting [22]. 

It should be acknowledged that the global extension of the singular field beyond the neighborhood of the North 
Pole is not unique. This ambiguity, however, is of no consequence, because a different extension 3>'^ that continues 

to respect the local expansion through order e 1 will be such that <&^\ — n = 0(e 2 ). Because the difference vanishes 
at the position of the particle, the value of the self-force, and the value of the four rcgularization parameters, will not 
be affected. 

D. Squared-distance function 

An important piece of as can be seen from Eq. (2.62), is s 2 , the squared distance between the points x and 
x. This is defined by Eq. (2.19), and its leading term in an expansion in powers of w a is obtained by substituting 
a a {x,x) = —g a pu)P , a truncated version of Eq. (3.19), into Eq. (2.19). The result is 

P 2 ■= (9a/3 + u a up)w a w 13 , (5.18) 

in which the metric and the velocity vector are evaluated at x = [to, ro, § , <j>o}- With u a := [E/f , f , 0, L/r 2 ] we 
obtain 

We re-express this result in terms of w 1 ^ := (j) — 4>' = + cA and choose c in order to eliminate the term proportional 
to Awf With [19] 

C:= (r a -2M)(r 2 +L 2 ) (5 ' 20) 
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we find that Eq. (5.19) becomes 

where we have also used the geodesic equation to eliminate Tq, the square of the radial velocity at x, in favor of E 2 . 

We define a "squared-distance function" p 2 by making the substitutions w B = —\/2Q sinf!) and w'^ = \/2Qcos/3 
into Eq. (5.21); these are truncated versions of Eq. (5.16) and (5.17), respectively. We find that this is given by 



riE 2 

P := (^ 

in which 



X :=l-fcsin 2 /3, k :=^—^ (5.23) 

contains the dependence of p 2 on /3. Recalling the definition of Q from Eq. (5.15), the squared-distance function can 
also be written as 

p 2 = 2(r 2 + L 2 ) X (S 2 + 1 - cos a), (5.24) 

with 

K 2 t 4 A 2 

& .= tLlo £±_ Cco^ 

• 2(r 2 + L 2 ) 2 (r - 2M) 2 X ' 



E. Calculation of s — continued 

At the end of Sec. V C wc had computed the singular field * (to, ro + A, a, (3) and expressed it in terms of sin /3, 

cos (3, Q := \/l — cos a, and A; the expression also involves world-line quantities such as ro, E, L, and ro which 
describe the state of motion of the particle at x. We interrupted this computation in Sec. V D to introduce the 
"squared-distance function" p 2 , which is the leading term in an expansion of s 2 := (g a/3 + u a u l3 )a a <r^ in powers of Q 
and A; this is itself a function of sin/3, Q and A, as well as world-line quantities. 

The squared-distance function is introduced to eliminate the dependence of on all even powers of Q, and to 

reduce all odd powers of Q to something linear in Q. The idea is to solve Eq. (5.22) for Q 2 and to substitute the 
resulting expression Q 2 (p 2 ,A 2 ) into our current representation of the singular field. We thus systematically replace 
each factor Q 2n in $^ by [Q 2 (p 2 , A 2 )]", and each factor Q 2n+1 by [Q 2 (p 2 , A 2 )] n Q. This yields a representation of 
the singular field which separates into a first set of terms that is independent of Q, and a second set of terms that is 
proportional to Q; each set contains a dependence on sin/3, cos/3, p, and A, as well as on the world- line quantities. 

Our final expression for $? J to, ro + A, a, (i) is too long to be displayed here. In fact, it is too long to be displayed 
anywhere, and we have taken measures to keep it hidden within the depths of our Maple worksheets. Its basic 
structure is as follows. The singular field admits an expansion in powers of e of the form 

= *(m),-2 + + $ (Uo + *?„),+! + 0(e 2 ), (5.26) 

in which $^ _ 2 [ s f order e -2 , $^ _ 1 of order e _1 , and so on. (Recall from Sec. II that e loosely measures the 
distance between x and x; we have that p and A are both of order e.) The terms that appear on the right-hand side 
of Eq. (5.26) possess the schematic form 

*(„),-2 - Af (A1 ),- 2 (A/p 3 )+0(Qcos/3/p 3 ), (5.27) 
*(„),-! = M (Al)i _ 1 (l/p) + 0(Qcos/3A/p 3 ) + 0(A 2 /p 3 ) + 0(Qcos/3A 3 /p 5 ) + 0(A 4 / P 5 ), (5.28) 
$ („),o = 0(Q cos /3/p) + 0(A/p) +0(Q cos /3A 2 /p 3 ) +0(A 3 /p 3 ) +0(g cos /3A 4 /P 5 ) 

+0(A 5 /P 5 ) +0(Q cos /3A 6 /p 7 ) + 0(A 7 / P 7 ), (5.29) 
= M (Al)!+1 p + 0(Qcos/3A/p)+0(A 2 /p) + 0(gcos / 3A 3 /p 3 )+0(A 4 /p 3 ) 

+0(gcos/3A 5 /p 5 ) + 0(A 6 /p 5 ) + 0(Q cos /3A 7 /p 7 ) + 0(A s /p 7 ) 

+0(Q cos (3A 9 /p 9 ) + O(A 10 /p 9 ). (5.30) 
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These equations display the dependence of each term on Q, cos /?, A, and p; the remaining dependence on (3 is 
contained entirely in the function \ '■— 1 — A: sin 2 (3 defined by Eq. (5.23). This dependence, as well as the dependence 
on the world-line quantities, is left implicit in Eqs. (5.27)-(5.30). The terms in these equations that involve the 
coefficients Mt^_ 2 , Mt^_i, and M^ +1 arc important: As we shall see, only these terms actually contribute to 
the regularization parameters. The coefficients depend on \ an d the world-line quantities. The remaining terms in 
Eqs. (5.27)-(5.30), those represented by the various symbols 0( ), are unimportant: They do not contribute to the 
regularization parameters. 

F. Final results: Regularization parameters 

The singular field of Eqs. (5.26)-(5.30) can now be substituted into Eqs. (5.10) and (5.11) to compute the multipolc 
coefficients The techniques reviewed in the Appendix provide us with efficient calculational rules. Equation 

(A. 34), for example, implies that all terms involving cos (3 in Eqs. (5.27)-(5.30) average to zero and do not contribute 
to the multipole coefficients. Equations (A.38)~(A.42), on the other hand, imply that all remaining 0( ) terms in 
Eqs. (5.27)-(5.30) vanish in the limit A — > 0; recall that this limiting procedure was introduced back in Eq. (5.2). 

The only surviving contributions come from the terms involving the coefficients M^_ 2 , M^- ) __ 1 , and M^ +1 . 
These are handled with the help of Eqs. (A.35)-(A.37), and this shows that does indeed take the form of 

Eq. (5.6). We remark that the neglected 0(e 2 ) terms in Eq. (5.26) are those which would produce the neglected (• • •) 
terms in Eq. (5.6); the (• • •) terms sum to zero because the 0(e 2 ) terms in the singular field vanish at the position of 
the particle. 

The regularization parameters are produced by the term involving the coefficient M^_ 2 in Eq. (5.27). Our 
results are listed back in Sec. I E, in Eqs. (1.30)-(1.32). Notice that in these equations, we changed our notation from 
ro to r(t), and from 0o to <p(t). The regularization parameters Bi^ arc produced by the term involving the coefficient 
M( M ) t _i in Eq. (5.28). Our results are listed back in Sec. 1 E, in Eqs. (1.33)-(1.37). The regularization parameters C( p ) 
would normally have originated from Eq. (5.29). Because, however, there are no surviving contributions from , 
we conclude that CV^) = 0. The regularization parameters D(n) are produced by the term involving the coefficient 
M( m))+1 in Eq. (5.30). Our results are listed back in Sec. I E, in Eqs. (1.41)-(1.45). 

The regularization parameters of Eqs. (1.30)-(1.45) are expressed in terms the (rescaled) elliptic functions S := 
f Jo" /2 (l-fcsin 2 V) 1/2 # = F{-±, i; 1; k) and X := f £ /2 (l-ksin 2 V)" 1/2 # = F{\, §; 1; k), where k := L 2 /(r 2 + 
L 2 ) was introduced in Eq. (5.23). As indicated, the elliptic integrals can also be expressed in terms of hypergeometric 
functions. These appear naturally in the course of averaging over the angles (3 — see Eq. (A. 28). In fact, to obtain 
our final expressions for the regularization parameters we have rationalized their dependence on the hypergeometric 
functions by using the recurrence relation of Eq. (A. 31). 

As a final remark we note that the regularization parameters depend on A in two distinct ways. First, Eqs. (1.30)- 
(1.32) show that the parameters Ai^ are proportional to sign(A) and therefore discontinuous across A = 0; this 
behavior accounts for the discontinuity across r = r of each mode <&^)i m of the retarded field. Second, the reg- 
ularization parameters contain an implicit dependence on A that is not shown in Eqs. (1.30)-(1.45). Indeed, the 
right-hand side of each equation should include terms that depend on A but vanish in the limit A — > 0; they are not 
displayed precisely because they do not survive the limiting procedure described in Eq. (5.2). 
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APPENDIX: MULTIPOLE COEFFICIENTS 

In this Appendix we collect results from the literature that are required for the computation of the regularization 
parameters in Sec. V. Mostly we rely on the techniques developed in the paper by Detweiler, Messaritaki, and Whiting 
[22]. 
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1. Decomposition of a scalar function in spherical harmonics 

Let f(9, 4>) be a scalar field on S 2 . We decompose it in spherical harmonics as 

f{6,<t>) = Y,flrnY lm {6,<t>). (A.l) 
Im 

The sum over the integer I extends from I = to I = oo, while the sum over the integer m ranges from m = — I to 
m = I. The field's spherical-harmonic modes are given by 

fim = J f(9,cj>)Y lm (8,cf>)dn, (A.2) 

where an overbar indicates complex conjugation, and dfl = sin 9 d9d(f> is an element of solid angle. 
We shall be interested in the value of / at the special point (6 = <f> = 0). This we express as 

/(§,0) = £/i. (A.3) 
i 

with 

l 

fi= £ .f Jm y' m (f,o). (A.4) 

m=— I 

The quantities /; associated with f(9, </>) play a fundamental role below; we shall refer to them as the multipole 
coefficients of the function /. 



2. Rotation of the angular coordinates 



A convenient way to calculate the multipole coefficients fi is to perform a rotation of the angular coordinates that 
maps the special point (§,0) to the North pole of the new coordinate system [17, 18]. The rotation from the old 
angles (6, <j>) to the new angles (a, (3) is defined by the equations 

sin 8 cos 4> = cos a, (A. 5) 

sin 9 sin (f> — sin a cos /?, (A. 6) 

cos 9 = sin a sin fi. (A. 7) 

It is easy to see that this does indeed map the point (9 — |, <f> = 0) to the point (a = 0, [3 =?), with /? undetermined. 
The special point is a singular point of the new coordinate system, and as we shall see, this property is a source of 
simplification in the computation of /; . 

It is well-known that a rotation of the angular coordinates changes Yi m (9, 4>) into a mixture of functions Yi m i (a, (3) 
with m' ranging from —I to I; the rotation mixes m but leaves I invariant. Because the quantity on the left-hand-side 
of Eq. (A.3) is a scalar, whose value is unchanged by the rotation, we may conclude that the rotation leaves the 
multipole coefficients /; invariant. We shall make use of this simple fact to find an efficient way to compute these 
quantities. 

The function /(a, (3) can be expanded in spherical harmonics Y lm (a, [3) as in Eq. (A.l), and such a decomposition 
implies that /(0, ?) = ^2 im fi m Y lrn (a = 0, (3 =?) = J2i V (2£ + l)/(47r)/zo, after evaluating the spherical harmonics. 
Comparing this with Eq. (A.3) yields 



/' = Y ~4~~ -fo[( a 'ft) decomposition]. (A. 8) 

As indicated in Eq. (A. 8), the multipole coefficient /; is proportional to the axisymmetric mode fio of the function 
f(a,[3). This is given by 



f^ = \j-^- J f(a,(3)Pi(cosa)dcosad/3, (A.9) 
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where Pj(cosa) is a Legendre polynomial. Equation (A. 8) therefore becomes 

ft = Hi±J: [ f( a ,/3)Pi(cosa)dcosadp. (A.10) 

47T 7 



If wc let 



1 f 1 

/,(/?) :=-(2/ + l) J J(a,P)Pt(cosa)dcosa, (A.ll) 
then Eq. (A. 10) can be expressed as 

fi = (MP)), (A.12) 



where 



1 



27i 



CM/3)> = ^ ^ //(/W (A.13) 

is the average of fi{(3) over all angles (3. Equations (A.ll) and (A.12) summarize the method by which the multipolc 
coefficients of the function f(0, (f) are computed. 

Equations (A.11)-(A.13) can be given an alternative interpretation that turns out to be useful for our purposes. 
Suppose that we are presented with a scalar function f(a,(3) and that we wish to represent its dependence on a in 
terms of an expansion in Legendre polynomials. We would write 

f(a,p) = J2MP)Pi(cosa), (A.14) 

and we would use the orthogonality properties of the Legendre polynomials to express /;(/?) as in Eq. (A.ll). The 
average of fi((3) over the angles (3, as defined by Eq. (A.13), would then give us the multipole coefficients ft. It is this 
interpretation that will be emphasized in the rest of this Appendix. 

We pause here to remark that the calculational method described above to compute the multipole coefficients /; 
takes advantage of the fact that the rotation of the coordinate system maps the special point (6 = ^ , <j> = 0) to the 
singular point (a = 0, (3 =?)• It is the singular nature of the new coordinates at the North pole which produces the 
equality of Eq. (A. 8). This method works well because f(9, (f>) is a scalar function of the angles: the rotation does not 
change the numerical value of the function. The method would not work as well for vectorial or tensorial functions of 
the angles, because the map to a singular point would alter the value of the functions by a singular factor. While this 
difficulty can be averted [17, 18], this introduces complications that can be avoided by choosing to deal with scalar 
functions only. It is this observation that has motivated us to work in terms of the scalars := instead of 

the vector $ Q := V Q $. 

3. Distance function p(a, (3) 

The multipole coefficients that are required in Sec. V are those associated with negative and positive powers of the 
"distance function" p(a,(3), defined by Eq. (5.24), 

p 2 :=2{rl+L 2 ) X {5 2 + l-cosa). (A.15) 

The symbols that appear in Eq. (A.15) are all introduced in Sec. V D. We have ro denoting the current radial position 
of the charged particle, and L is the particle's (conserved) angular momentum per unit mass. In addition, 

a2 - 2 , T^JZ 0,^2 — (A.16) 



■ 2{r 2 + L 2 ) 2 (r - 2M) 2 X ' 
in which E denotes the particle's (conserved) energy per unit mass, and 

A:=r-r , (A.17) 
is the radial component of the displacement vector w a = x a — x a . And finally, 

X^l-fcsin 2 /?, fc:=^L^ (A.18) 

contains the dependence of p on the angle {3. 
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4. Decomposition of (S 2 + 1 — cos a) n in Legendre polynomials 

The calculation of the regularization parameters described in Sec. V requires multipole coefficients for various powers 
of p, in a context in which A, and therefore S, is small. The relevant powers of p are p~ 2 ™ -1 with n = — 1, 0, 1, 2, and 
so on. As was explained in the paragraph surrounding Eq. (A. 14), the starting point of a computation of multipole 
coefficients is the decomposition of those selected powers of p in Legendre polynomials. The most important piece of 
p for these decompositions is the factor (S 2 + 1 — cos a) 1 / 2 , and we therefore need expansions of the corresponding 
powers of this quantity in Legendre polynomials. 

Quoting from Detweiler, Messaritaki, and Whiting [22], we have 

(<5 2 + 1 - cos a)-™" 1 / 2 = J2 A ?( S ) p i( cosa )' ( A - 19 ) 

i 

in terms of coefficients Af(6) that can be expanded in powers of S. They obey the recurrence relation 

• 4 " +1 " "(STijisk- < A ' 20 > 

and special values are given by 

A° t = V2-(2l + 1)6 + {21 + }) 2 S 2 - h(l + 1)(2Z + l)5 3 + 0(<5 4 ), (A.21) 

2V 2 3 

A} = ^- { 21±}L + , ( , + im + l)S + 0(S>), (A.22) 

^.^.ffiifil!^, (A.23) 
A 3 = ^_i™±i) +0(rI)> (A , 4) 

A . = n^ji ± m ±21+0{5 ^ (A , 5) 

By induction from Eq. (A. 20) we infer that 

\l + 0(5 2 )}, n>2. (A.26) 



A- 21 + 1 



2n-l 



{2n-l)5 : 

We shall also need a decomposition for (<5 2 + 1 — cos a) 1 / 2 . This is given by Eq. (A. 19) with n = —1, and in this 
case the expansion coefficients are [22] 



2\/2 
{21 - 1)(2Z + 3) 



A- 1 = - 7 —JX±—^ + 0(5 2 ). (A.27) 



The results displayed in this subsection can be used in concert with Eq. (A. 15) to calculate the Legendre coefficients 
jp-2n-i]j^ associated with the functions p^ 2n ^ 1 (a, [}). The dependence on [3 is contained in the factors of \ that 
are hidden in the definitions of 6 and p; refer back to Eqs. (A. 16) and (A. 18). 



5. Averaging over f3 

The next step in the calculation of the multipole coefficients (p~ 2n ~ 1 )i is to carry out the averaging over the angles 
(3, as defined by Eq. (A. 13). Because the dependence on [3 is contained in x, what we need are expressions for (x~ p ), 
where p is some (positive or negative) number. 

Again quoting from Detweiler, Messaritaki, and Whiting [22], we have 

{x- p )=F P ~F{p,\;l;k), (A.28) 

where F(a, b; c; x) is the hypergeometric function, and where k is the constant defined by Eq. (A. 18). Special cases of 
Eq. (A.28) are 
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and 

(X) = 1 - \k. (A.30) 

The hypergeometric functions F p are linked by a recurrence relation displayed in Eq. (15.2.10) of Ref. [41]. When 
specialized to our specific situation, this equation reads 

Using this identity we find that the F p 's that appear in our expressions for (p~ 2n_1 )z below — see Eqs. (A. 36) and 
(A. 37) — can all be expressed in terms of F 1 / 2 and F_x/ 2 . These, in turn, can be expressed in terms of complete 
elliptic integrals: According to Eqs. (17.3.1), (17.3.9), and (17.3.10) of Ref. [41], we have 

F 1/2 = X := -K(k) := - / ' (1 - k sin 2 ipy 1 ' 2 dip (A.32) 

and 

F_i /2 = S := -E(k) := - [ ' (1 - k sin 2 ip) 1 ' ' 2 dip. (A.33) 

7T IT J Q 

Our results in Sec. V E are expressed in terms of the elliptic integrals. 

We conclude this subsection with the remark that any function g of x is necessarily periodic with period tt in the 
interval < (3 < 2n; the function's behavior in the interval < (3 < n is replicated in the interval tt < (3 < 2n. 
Furthermore, any such function is necessarily symmetric (even) about (3 = ir/2 in the first interval, and about (3 = 3ir/2 
in the second interval. These properties are sufficient to infer that 

(<?(x)sin/3) = (<Kx)cos/3) = (g( X ) sin /3 cos/?) = (A.34) 

for any function g(x)- 



6. Final results: Multipole coefficients (p 



The multipole coefficients of the functions p~ 2n_1 (a, (3) are denoted (p^ 2 " -1 );. They are calculated by implementing 
the procedure described around Eq. (A. 14). Most of the required computations were already performed in Sees. 4 
and 5 of this Appendix; here we simply collect the results and put it all together. 

By combining Eqs. (A.15), (A.16), (A.19), (A.21)-(A.27), as well as Eqs. (A.28) and (A.29), we arrive at the 
following listing of multipole coefficients: 



r - 2M 



Erl 



sign(A) + 0(A), 



p+l/2 



L 2 



+ 0(A), 



(x- p P )i 



" +1) 



(i-m 



+ 0(A 2 



We also record the scaling relations 



which follow from the same set of equations. 



A(x~ 


V P 




= 0(A), 


A 2 (X- 


- p p- 




= 0(A), 


A 4 (x~ 


- p p- 


" 5 ) ; 


= 0(A), 


A 7 (x 


- p P 


- 7 ) ; 


= 0(A 2 ), 


A 10 (x 


P P 


9 ) ; 


= 0(A 3 ), 



(A.35) 
(A.36) 

(A.37) 



(A.38) 
(A.39) 
(A.40) 
(A.41) 
(A.42) 
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